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THE TRUNCATED POISSON 
F. N. Davin anv N. L. JoHNSON 
University College, London 


1. Distributions often arise in practice which are of the Poisson 
type but in which the zero group is unobserved. Examples of such 
cases may be found in the distribution of the number of persons per 
house suffering from an infectious disease such as measles, or in the 
distribution of the number of accidents per worker in a factory. It 
was the latter distribution which brought the problem to our notice, 
for while it was a simple matter to count over a given period of time the 
number of workers sustaining one, two, or more accidents, the number 
of persons who did not have an accident could not be enumerated 
owing to the factory population fluctuating during that period. We 
shall discuss here problems of discrimination and estimation connected 
with distributions of this type. 


2. Given a quantal distribution with no zero group observed, 
before the number in this zero group is estimated, it is necessary to 
decide about the functional form of the distribution. When the zero 
group is known the simplest test to use to discriminate between a simple 
and a compound Poisson hypothesis is the variance over mean ratio. 
If the observations be x, , x. , -*- , ty , then we calculate the index of 
dispersion, 


t=1 
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and refer to x° tables with (NV — 1) degrees of freedom. Kathirga- 
metamby (1952) has shown that this approximation leads to actual 
significance levels between 4.96% and 5.1% if a nominal signifieance 
level of 57 is used, even if N is as small as 20, provided \, the parameter 
(mean) of the Poisson distribution, is not smaller than unity. This 
procedure seems likely to be satisfactory for most cases encountered 
in practice. 

It did not prove possible to devise a similar type of test for the 
truncated distribution which would be independent of the unknown 
value of \ for the complete Poisson distribution. We resort therefore 
to the procedure described below. 


3. Assume that we have N observations in the discrete distribution 
truncated by omission of the zero class, and that the unknown number 
in this missing class is nm) giving a total of (NV + mn.) observations in 
all. Letting primes refer to values obtained from the complete sample 
the index of dispersion for the completed distribution will be 


= = 


(N + 


(N + 2 — Ne. 


Ni 


The significance of J’ is judged by reference to the x’ tables with (V + 
N — 1) degrees of freedom. Since usually N + n,. — 1 will be large 
this means approximately that the hypothesis of a Poisson distribution 
will be rejected when 


2 - 
Ni N + N + 
where k is a positive number depending on the level of significance of 
the test. From the inequality we see that 
(a) there will always be a number N, such that the Poisson hypoth- 
esis would be rejected if no > No , 
(b) if the Poisson hypothesis would be rejected with n, = 0 it 
will be rejected whatever the value of np . 
In this latter case it seems a safe procedure to regard departure from 
the Poisson distribution as established. It is easy to see that this 
condition is satisfied if the index of dispersion is applied to the actual 
data, ignoring the effect of truncation, and a significant result is obtained. 
It must be noted that such a procedure will, in general, lead to the 
actual probability of rejections being less than the nominal probability. 
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4. Numerical Examples. We now give two numerical examples of 
the point mentioned in §2._ We take complete data and apply the Index 
of Dispersion test both before and after truncation of the zero group. 


TABLE 1. DEFECTIVE TEETH IN 11 YEAR OLD BOYS 


Number of 
Tecth 

Number 
of Boys | 61 | 47 | 43 | 35 | 28 | 15 | 20| 2] 2 1 | 265 


| 
@ 


10 | 11 | 12 |Total 


For the complete data the mean is 2.574 and the Index of Dispersion 596. 
This latter is equivalent to a unit normal deviate of 14.4 and is obviously 
highly significant. Thus the data will not be adequately described by a 
Poisson distribution and we should look for a suitable compound Poisson. 
We now leave out the zero group and consider 204 observations at points 
1, 2,3, --- . The mean is 3.343 and the Index of Dispersion 302, giving 
a unit normal deviate of 4.95. Again we note that this is highly signi- 
cant although much less so than that obtained from the complete data. 


TABLE 2. TICKS ON SHEEP 


Number of Ticks 6| 8| 91 10 jTotal 


Number of Sheep 


on 


R. A. Fisher quotes the above data concerning the number of ticks found 
on sheep. The Index of Dispersion on the total data is 107.5 which on 
reference to x” tables, f = 59, is significant at the 0.5% level. Omitting 
the seven observations at zero the Index of Dispersion on the data 
thus truncated is 72.16 and the degrees of freedom 52. This is signifi- 
cant at the 5% level. 

The method of approach so far employed avoids discussion of one 
very important difficulty. Are we to regard our (truncated) sample as 
one of a sequence of samples of equal size (i.e. N constant) or as a result 
of truncation of one of a sequence of complete samples of the same 
size (NV + n, constant)? 

In the following analysis we shall take the first of these viewpoints 
which leads to rather more straightforward sampling theory than the 
second, but it must be realised that our analysis may well be inapplicable 
in certain problems. 
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5. We now suppose tests have shown that the Poisson hypothesis 
is a reasonable one to entertain with regard to the truncated distribution. 
There are various points of estimation which we may now wish to 
consider but they all reduce to the basic problem of estimating the 
unknown parameter of the Poisson series which we shall call \. We 
write small letters for the sample moments about the origin and capital 
letters for the corresponding population moments. It is obvious that 


A= M’ 
and accordingly the simplest estimate that we may calculate is 
, 
m, 


Approximations to the moments of log (\ + 1) are easily found, at 
least to order 1/N’ , where N is the number in the truncated sample. 
For the first moment, for example we write 


log (X + 1) 


i! (- ar) + 


and we average over all samples in which there are N observations 
after truncation. We have then, to order 1/N, 


2N 


or, in terms of the unknown parameter }, 


log + 1) = log +1) 
Fiog(i+1) = NO + 1? 


& log (A + 1) = log 


it is easy to show from this that 


N 


Thus, provided N be of a reasonable size the estimate should be good 
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enough. We illustrate this point by means of table 3, giving values 
of A and B such that, to order N~ . 


A 
& log (A + 1) = log (A+ 1) 
2 B 
Flog(i+1) = 
TABLE 3. 
N A B log (A + 1) 
0.5 0.48091 0.43719 0.17609 
1.0 0.55311 0.47409 0.30103 
1.5 0.52827 0.43505 0.39794 
2.0 0.48037 0.38430 0.47712 


6. The maximum likelihood estimate of \ in the truncated case is 
difficult to find but we may obtain an expression for the variance of this 
estimate. The rth term of the truncated distribution is 


whence the joint probability function is 
= \** constant. 
(1 


where Z is the truncated sample mean and N, as before, is the number 


in the truncated sample. Equating the differential of L = log P to 
zero we have 


Ne Ni 
i.e. ey 


from which it does not seem possible to obtain an explicit expression ie 
‘for \. For the variance of \ we have 


1 1 Al 


From this latter we may make a comparison of the efficiency of the 
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crude moment method with that of the maximum likelihood method. 
Thus 


Var. (A) Max. Lik... 1 
Var.(X) Moments *A+2 1 —a/e — 1) 


this ratio varies with \ and has a minimum value between \ equal to 
to 2 and A equal to 4. 


TABLE 4. RATIO OF VARIANCES OF MOMENTS AND MAXIMUM 
LIKELIHOOD ESTIMATES 


r 
Ratio 


NI 


1.0 1 
7 0.80 0. 


The ratio will never be less than 70% and with increasing ) will tend 
to 100%. It would appear therefore that while the method of maximum 
likelihood will give estimates which have smaller variance than the 
method of moments, the ease with which the latter estimate is calculated 
may be thought to outweigh the advantages of the former, particulafly 
since the efficiency ratio never drops below 0.70. 

On the other hand the maximum likelihood solution will be the 
less troublesome if appropriate tables are available. Table 5, giving 
X (Maximum Likelihood) corresponding to @ = 1.1 — (0.1) — 3.5 will 
be of some help. For # > 3.5, = 2(1 — e*) is a good approximation 
(error les. than 0.1%). 


TABLE 5. 

z i z 
1.1 0.194 2.1 1.726 3.1 2.935 
1.2 0.376 2.2 1.856 3.2 3.048 
1.3 0.550 2.3 1.984 3.3 3.160 
1.4 0.715 2.4 2.109 3.4 3.271 
1.5 0.874 2.5 2.232 © 3.5 3.381 
1.6 1.027 2.6 2.353 

1.7 1.175 2.7 2.472 

1.8 1.318 2.8 2.589 

1.9 1.458 2.9 2.705 

2.0 1.594 3.0 2.821 


7. We turn now to the case when the hypothesis of a Poisson dis- 
tribution is rejected, and the physical conditions of the problem suggest 
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that a compound Poisson might be a suitable alternative. For illustra- 
tive purposes we shall choose the negative binomial 


where p= and s are positive quantities and gq = 1 + p, as our compound 
Poisson. We shall show that to use the method of moments for the 
estimation of p and s leads to very inefficient results by any standards 
but especially when compared with the maximum likelihood estimate 
variance. We shall derive this latter estimate first. The probability 
of observing r successes may be written as 


P, = -q’) r=1,2,-- 


Writing 7 for the sum of the observed values, we have the likelihood 
function 


+ T log p — (Ns + T) log q — N log (1 — g’) 


from which we derive the two equations which are, theoretically, soluble 
for the estimates of s and p. Thus 


s+r-1 
op @ 
PQ. — 


There appears to be no way of obtaining estimates of p and s except by 
iteration. The variance of the estimate of p, say p, is approximately 
given by 


Var. (p) = 


af ast] ® op” 3s op 
Further differentiation of the equations above shows that 


_ Na “(log 9)" 
as" 
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2 

aL = ie = 1) +— Ns 

We notice that d°L/ds dp does not depend upon r or upon T and so it is 
unaltered by the operation of taking expectations. Again since 


&(T) = N&(r) = Nsp/(1 — 


it follows that 


(54) = +@~- Dg") . 


Rather more trouble is involved in the calculation of the expectation of 
0°L/ds° and the best that can be done with it appears to be to express 
the first term as a series. Thus 


_ Nq “(log q)’ 
1 


N —s ] 2 

The arithmetic involved in the evaluation of these various quantities 


is not difficult and the variance of the estimate of p proves surprisingly 
small. Thus we have 


1 


TABLE 6. VALUES OF VARIANCE OF ESTIMATES OF p MULTIPLIED 
BY THE TRUNCATED SAMPLE SIZE 


& 

om 


N Var(p) 1.178 1.680 


8. Estimates in terms of the sample moments are often to be pre- 
ferred over maximum likelihood estimates in that they lead to simpler 
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equations to solve, and ease of calculation will often outweigh loss of 
efficiency. In the maximum likelihood equations we see that functions 
appear involving only the sum of the observed values. For moment 
estimates it is necessary to use the first three moments o! the truncated 
sample. It is this introduction of the third sample moment which 
makes the method of moments particularly inefficient in this connec- 
tion. The [th moment of the truncated negative binomial may be 
written as 


— (s — 1)! 


This leads one to suggest that we should take for an estimate of p, 


M; = 


We expand about the point M{ , Mj , M$ and, after taking expectations, 
obtain to order 1/N, 


_ Mi 2Mi + MiM; — 
(M; — 


Substitution for the M’s in terms of s and p gives 


~ 


+ p(s’ + 12s + 19) + 4(s + 2)] 


The bias is large and we notice that as p, or s, becomes larger this 
bias in the estimate of p, for a fixed truncated sample size, becomes 
worse. Again 


_M Mi — M3 M3 Mi — M 1 
+ an — my 7M {its uy t 


J Mi Mi, | 1 
+ uy t 
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Substitution for the crude moments in terms of s and p does not seem 
worth while since a simple example will serve to show how very poor the 
variance of this estimate is. Thus if we put s = 1 = p we obtain 


N Var (p) = 69 


i.e. a standard error which is some eight times that obtained by the 
maximum likelihood method. In fairness to the method of moments 
we should point out that the largeness of this standard error is almost 
entirely due to the introduction of the third sample moment, but again 
the only point in taking the method of moments is in order to use as 
many moments as will give an early solution. We would recommend 
therefore a process of iteration with the maximum likelihood equations 
rather than a moment solution for this particular problem. 

Such an iterative procedure is outlined below. The equations for 
p(= q — 1) and ¢ are 


+ 1) — NWG) = (N log — 
T = — 
where (zx) is the digamma function (Davis 1933)). 
Hence 
(log = 3{ Ws +r) — 
For a given value, %, , of 8, we can determine a corresponding value, 


Dr , of p, using this equation. A new estimate of s, s, , can then be 
calculated as 


8 = T(1 — . 
% in its turn can be used to find a new estimate, Do , for D, and soon. 


Table 7 is intended to help in the solution (for p) of the equation 
(log g)/p = y. 
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TABLE 7. 

y y p 
0.40 4.047 0.75 0.734 
0.45 3.177 0.80 0.538 
0.50 2.513 0.85 0.373 
0.55 1.994 0.90 0.230 
0.60 1.579 0.95 0.107 
0.65 1.242 1.00 0.000 

0.70 0.965 
REFERENCES 


H. T. Davis (1933). Tables of the Higher Mathematical Functions, Vol. I. Principia 
Press, Bloomington, Indiana. 
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INVERSE, MULTIPLE AND SEQUENTIAL 
SAMPLE CENSUSES’ 


Dove.as G. CHAPMAN 


University of Washington 


psa enumeration of populations by sampling methods has only 
recently come into widespread use, though the idea dates back at 
least to 1783 when it was proposed by Laplace [9]. Approximately a 
century later, Petersen [13] used the tag and sample method to enumerate 
a plaice population. The basic procedure consists of tagging or marking 
a number of individuals within the population and subsequently samp- 
ling, at random, the whole population. Estimates of the population 
are based on the random number of tag recoveries in the sample; the 
number of individuals marked and the sample size are known parameters. 
A study of such estimates and of tests for the population size has been 
made by the author {3}, [4]. 

While this simple model suffices for many purposes it is apparent 
that more complex models will be often necessary or desirable. In 
many cases the tagging and sampling is carried out in several stages. 
Such a procedure is referred to, here, as a multiple sample census. 
Various point and interval estimation formulae based on this procedure 
have been given by Schnabel [17], Schumacher and Eschmeyer [18], 
De Lury [7], [7a], Schaeffer [16] and the author. This type of procedure 
lends itself to a sequential procedure. Tests and estimation formulae 
based on such multiple and sequential procedures are formulated in this 
paper. 

An even simpler modification of the single step census is to “invert” 
the sampling process i.e. to fix the number of tagged individuals to be 
recovered by sampling, rather than fixing the sample size. Some 
formulae using this idea were given recently by Bailey [1]. 


1Research partially supported by the Office of Naval Research. 
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INVERSE SAMPLE CENSUSES 


The following notation will be used throughout: 
N: the population size being studied 
t,: the number of marked individuals in the population at the time 
the 7-th sample is taken 
n,: the number of individuals taken in the 7-th sample 
s,: the number of tagged individuals recovered in the 7-th sample 
s;,: the number of individuals recovered in the i-th sample that 
had been tagged at the j-th tagging. 
For convenience we also define n» = ¢; and s, = 0. The subscript: 
i may be taken to run from 1 tok. Where k is 1, i.e. the sample census 
is conducted in one stage, the subscripts will be omitted for simplicity. 
If the inverted sampling plan is used, the s, are fixed parameters, 
the n, are random variables. For the single sample case n is a negative 
binomial or negative hypergeometric random variable according to 
whether sampling is with or without replacement. 
If sampling is with replacement, Pr(n), the probability of having 
to sample n individuals to obtain s marked ones, is given by 


and E(n) = sN/t. 

Hence nt/s is an unbiased estimate of N, which, as the author has 
pointed out [4], is not true when n is fixed and s the random variable. 
The variance of nt/s is (V?—Né)s"' and an unbiased estimate of this 
variance is 


2 _ nt(n — 8) 
(2) + 1) 


As noted in Appendix I 


7 nt/s — N 
VN°*/s — Nt/s 


(3) 


is approximately distributed according to N(0,1) for large s. This fact 
may be used to set up confidence intervals and tests for N. 

For s and ¢/N both small, as will be more frequently the case, the 
second limiting distribution given in Appendix I will be more useful. 
For this range of the parameters 2nt/N has approximately a x’ distri- 
bution with 2s degrees of freedom. This is equivalent to the Poisson 
approximation to the binomial and is related to the results of Sandelius 


(1) Pr(n) 
: 
* 
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[15] concerning inverse sampling with random variables distributed 
according to a Poisson distribution. 
If x3. (€) denotes the e-quartile of the x’ distribution with 2s degrees 
of freedom i.e. 
S = € 
then (1 — e) confidence limits for N are given by 


(4) 


The problem of improving on these confidence limits in a procedure 
similar to that treated in [3] will be considered elsewhere. 
If the sampling occurs without replacement 
(n — — — n)! 


and 


_ (N + Ds 
6) = 
so that 
(7) 


is an unbiased estimate of N. In contrast to the direct sample census 
estimate, this unbiasedness does not depend on the parameters s and t. 
Furthermore 


2 _ (N+ 1(N N 


This approximate formula (which exaggerates the actual variance) 
may be useful in the determination of the choice of s: by an appropriate 
choice of 8, o@/N can be fixed at any desired level. 

For testing purposes we note that N is approximately normal with 
mean N and variance given by (8) for large s. A model similar to that 
used by David [6] can be constructed to prove that as s tends to infinity 
n is asymptotically normally distributed. 

On the average the inverse sampling procedure is better than the 
direct sampling procedure. For if n’, s’, denote the fixed sample size 
and the random number of tags recovered in the direct procedure, and 
if ¢ is chosen equal to n’t/N, then 
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while 
¥ (n’ + 1) + 1) 


is the almost unbiased estimate in the direct sampling ease. Hence a : 
more efficient estimate is obtained with less average effort. 

On the other hand if the experimenter knows absolutely nothing 
about the possible population size, then by an improper choice of ¢ and 
s, H(n) may be extremely large. Moreover 


2 aN + IN (*) 
(1+ 4 2) 


is very large; this may he regarded as an undesirable feature of the 
procedure. 

These difficulties may be partly overcome by a modification of the 
inverse sampling plan as follows: the number of untagged individuals 
to be recaptured is predetermined, rather than the number of tagged 
individuals. In other words, n — s is chosen in advance of sampling: 
s and n are now both random variables, though completely dependent. 
For convenience write n — s = 6. 

(9) in — — OUN — a)! 


(n — ot —n + 6 — — — IN! 
is derived in the usual manner. However the obvious estimate again 
is no longer strictly unbiased. For 


(10) — §+ i) = i)= 
(n— 8) (UN —1—1+1—8) (N+ 


(N+ — 1 


This result. is analogous to that obtained by the author in [4]; by 
the same argument and formulae given there (pp. 145-146) it follows 
that the second term in the parentheses is negligible provided 16‘N > 
log N. Vor such values of 4, t, N the estiniute 


n(t 1) 
N = 
(11) 1 


: 
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has bias less than 1 in absolute value. 


To determine the variance of N, (s + 1)°* is expressed in an inverse 
factorial series as in [4] so that 


(V+ 1? = + 1) + + 2) 


2 +h) 
(s + 1)(s + 2)(s + at (s + 1)(s + 2)(s + Bs + Ht 


To evaluate the expectation of this latter series let n,;, = n(n — 1)- 
(n — 2) --- (n — 7+ 1) and observe that for any i <7 (j —i = ©) 


— Den) (VN — 0), 
(12) (s+ Den (I — m) 


ITere n, stands for j terms of a form similar to the one in the brackets 
of formula (10). This formula is derived by a direct summation exactly 
parallel to that used in obtaining (10). 

If 7, is neglected and the remaining terms on the right hand side 
of (12) are written q,; then 


(13) oF E(XY N? = (i+ + 23 + 2424 + ---) 


— (qe + Gis + 2qu — 2N —2— 
Table 1 gives a tabulation of o¢,,N for various representative values of 
N, t and 6, obtained by means of this formula. 


The average sample size required by this procedure can be determined 
simply from (9). In fact again by direct summation 


+9), 
(14) En+i = 
so that 

(14) E(n) = 
and 


—t4+2) 

Kither of these formulae ean be obtained from (6) and (8) by an 
interchange of N -—- ¢ and ¢ and by replacing s by 6. It is seen im- 
mediately that the tremendous variation possible in the earlier inverse 
sampling model is now eliminated. 


(16) .= 


| 
koe 
qs 
1 
Ny 
cer 
=¥ 
tothe 


SAMPLE CENSUSES 291 


Since 6 will usually be reasonably large the most appropriate ap- 
proximation to use for testing purposes is the normal distribution. In 
particular, it is desirable to work with n using formula (15) and (16). 
Writing ¢/(N + 1) = p’ and using a trivial modification of the approxi- 
mation of (16) 


is approximately N(0,1) and confidence limits for p’ (and hence N) are 
obtained from the quadratic equation 


(17) — p’) — = — p’). 
The approximation may be improved slightly if the exact formula for 
a, is used but this involves solving a higher degree equation. 


DIRECT MULTIPLE SAMPLE CENSUSES 


In the most usual type of direct multiple sample census, at each stage 
a group of individuals are drawn from the population. Those that are 
tagged are noted; those not tagged are tagged and then the whole group 
is returned to the population. Such samples may bg of any size. 

For this situation the appropriate model is characterized by the 
following characteristics: 


(a) =n, — 8; i=1,2,---,k 
(b) (= mo), , M2, , are parameters 
(ec) to, ts, te & are random variables 


It may be pointed out that previous authors have generally regarded 
the t; (¢ = 2,3, --- , k) as parameters and it is true that at any stage 
i of sampling ¢, is known and hence conditionally it is a parameter. 
This fact is the basis of formulae (23) and (24) below. On the other 
hand, in setting up the “a priori” probability model and in performing 
the overall preliminary planning one must be cognizant of the fact that 
at the beginning of the experiment (, , 4, , --- , & cannot be exactly 
determined or predicted. 

The probability model is given by the equation 


— 
(18) 8) = Il 


i=l (") 


— | (V — n)! 
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The remaining terms telescope because of (a). From formula (18) 
it is apparent that }°_, s, is a sufficient statistic for N. Moreover the 


maximum likelihood estimate of NV is easily derived from (18) [see 
Appendix II]. It is the solution of the equation 


In view of the fact that Vis much larger than any 7, a first approx- 
mation to the maximum likelihood estimate is 


Mat; 
A 


This formula is obtained by ignoring terms involving 1/N° or higher 
powers of 1, N in (19). If the cubic terms are retained the approximate 
equation to be solved for the maximum likelihood estimate is 


k f \ 3 

The larger root of this equation (which will be denoted V,) is the 
desired one. As shown in Appendix I, under certain conditions that will 
be frequently met with in practice, the maximum likelihood estimate 
will lie between N, and N, and furthermore the difference VN, — N 
will be relatively small compared to V. 

It is interesting to note that the same estimate is obtained by the 
method of moments. For, if }>)., s, is denoted by S, , then 


¢ 


(22) E(S,) I] (1 


i 


kel 


where Rk; = sum of products of , n, taken 7 at a time 
with all subscripts diferent. bither of these results is most casily 
obtained by induction. 

The distribution of SS, is difficult to obtain in any simple uscable 
form, so that it is difficult to evaluate the smal! sample properties of 
these estimates. It may be noted that these estimates are not strictly 
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comparable with those of Sebumacher and Eschniever, Schnabel and the 
author's own, referred to in the introduction.  Phese latter estimates, 
which are generally of a xo type, are derived on the assumption (im- 
plicitly or explicitly) that the number of tags is a fixed parameter, not 
a random variable. 

It is possible to write down simply an unbiased estimate which is 
vilid in both of these models (i.e. whether the ¢, swe random variables 
or parameters) viz 


be +t ED | 


provided the s, (or more precisely n, (,N) are sufficiently large [roughly, 
near 10 or greater in size]. For 


DEIN] = 


under the restrictions noted (those given in [4] p. 115). 

As a consequence of a theorem of Blackwell [2], however it can be 
said that an unbiased estimate for N based on S, has a smaller variance 
than N*, within the model in which the ¢; are cumulative and ¢;,, — 
=n; — 8. 

Where the ¢; are not random variables but fixed parameters, V* ts 
a desirable estimate (subject to the restrictions n,/,,N not too small) 
and furthermore the variance of N* may be computed from formula 
33 of [4] and the fact that 


+1 
li will frequently happen that the n, ¢,,N are too small to permit 
N* being an unbiased estimate. For this situation, where the ¢, are 
fixed, a slight modification of an estimate given by Schnabel [17] will 
be most useful viz: 


+ 


This is based on the fact that each s, has approximately a Poisson 
distribution with parameter hence has approxi- 


(25) | Np = 


| 
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mately a Poisson distribution with parameter 


E(N») = — © *) is easily derived. Furthermore will be 
certainly negligible if ae . ut; > N-log N. This is a much lesser 
restriction than that the same inequality hold for each n,¢; . 

Neglecting terms of the form \’e* (j = 0, 1, 2, 3) an approximate 
formula for the standard deviation of NV is derived, viz., 


2,6 
(26) = (1 +> + 


i=1 


The derivation is exactly parallel to the derivation of (33) in [4], through 
the use of an inverse factorial series. 

Many sample censuses will not fall strictly into either of these 
models: some of the unmarked individuals taken in a sample will be 
returned as marked members, but not all. <A full treatment of this 
case is clearly complicated. Where the n,f;/N are small, the estimate 
N, is probably quite satisfactory. 

The Poisson approximation will also be useful as a basis of con- 
structing tests and confidence intervals for N. The confidence limits 
given in [3] may clearly be used (with né replaced by }oi., njt;) if 
the ¢, are fixed parameters. Within the model primarily considered in 
this section this will probably be reasonably satisfactory also. 

If the experimenter observes s;; the number of tags recovered in 
sample 7, that were placed in the j-th tagging, rather than simply 
s.(s, = >>}-1 8,,), more information is apparently available. Actually 
the knowledge of s,; is of no value in estimating N, or in constructing 
tests for N. For the joint probability of the s;; is 


~(n; —8;\(N — 


ial 
(| 


(n, — s,)! 


From the form of this probability distribution it is seen that 
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is still a sufficient statistic for N. 

However the s,, may be used to test the assumptions that the 
sampling is random (i.e., the probability of including an animal in the 
sample is independent of its being tagged) or that the population is 
constant. 

Under the hypothesis tested, that NV is constant, given 47 = 1, 2, 3, 
-++ , ¢) the s,; have a multihypergeometric distribution, or neglecting 
the sampling without replacement, a multinomial distribution. How- 
ever the usual x* test does not follow immediately for the estimate of N 
will be determined from the random variables s,; which are not inde- 
pendent for different z. The author, however, conjectures that in view 
of the nature of the x’ test this is asymptotically negligible. Much 
more serious is the fact that in many cases at least, /'(s;;) are too small 
for a reasonable application of the test. 

As an alternative procedure a non-parametric test is suggested. 
An array of the form 


nt; = 


may be formed, in which each element is a random variable with ex- 
pectation 1/N. 

Within each row the random variables are independent. Between 
rows that is not the case: for the s,, , for fixed 7, are observations from 
the same sample. However the correlation between any two s,; , say 

lia tig N — 


“NNN-1 
which will be usually negligible. Moreover if [;(, — s.)]N is small 
compared to x; the conditional probability of s;, given s.3 is almost 
equal to the unconditional probability of s,, . 
In view of this the sign test suggested by Moore and Wallis [12] to 
test for randomness in a sequence of independent observations from 


| 
k 
| 
: 
bea 
3 


206 BDIOMEPRICS, DECEMBER 1952 


2 common distribution may be appropriate. The test is based on the 
statistic D, the number of negative signs in the sequence of successive 
differences of observations. Moore and Wallis tabulated the probability 
distribution of D for small values of », the mumber of observations. 
They conjectured and Mann [11] subsequentiv proved rigorously, that 
D is asymptotically normally distributed. In application of the test 
since 
n— 1 n+ i 


= - and 


n—- 
Vn ti 
is taken to be distributed according to V(O,1) for xn > 12. If the array 
(28) is considered us a single sequence of observations the n of the 
Moore and Wallis test is here equal to {&(/e + 1)] 2. 

In many cases the alternatives to randomness are essentially one- 
sided. For example, some that may be considered are: 

(a) the tagged individuals die off more rapidly or disappear so as 
not to be available for sampling 

(b) the tagged individuals disperse from the tagging location slowly 
and are more likely to be recaptured in the samples taken soon after the 
tagging rather than later 

(c) there is a marked change in the composition of the population 
due to natural processes or to migration. 

If any of these alternatives is true, the random variables in each 
row of the array will tend to increase. In this case a test based on 
the whole array as a single sequence has the following defect: if the 
alternatives are true, in each row the probability of a negative difference 
is greater than 3, but the probability of a negative difierence between 
the last element of any row and the first of the next row will be much 
less than }. 

To avoid this it is necessary to consider each row separately i.e. the 
array (28) may be considered as / sequences of observations decreasing 
in length from & to 1. A test of randomness may be made using the sum 
of the number of negative differences in these k sequences (actually 
k — 1 since no difference is obtainable from the last row). 

Let 


A= 2,4 --- + 


where D, = number of negative differences in row i. 
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Then 


(X) = — 1) 


_ (k§ + 4k — 1) 
4 4 


and oy 


The asymptotic normality of XY is immediate using the fact that D,_, 
itself tends to be normally distributed as / tends to infinity while the 
initial terms of the sum (29) are asymptotically negligible as k becomes 
large. 

Table II gives a partial tabulation of the distribution of X for 
k = 4, 5, 6, 7. 


INVERSE MULTIPLE SAMPLE CENSUSES 


As in the single stage eensus two models may be considered. In 
ihe first. of these, sampling and tagging are carried out in exactly the 
sume manner as in the direct multiple sample census except that the 
number of tagged individuals at each stage is predetermined rather than 
the sample size. 

In this mode] 

(a) —-t; =n, — 8; +=1,2,---,k 

(b) t, (= no), , , , & are parameters 

tp, ts, m are random variables, and 


Prin, , My) 


Il (n; — — — n,)! 
(8; — — 8) — s)INYUN — n, — 


There is now no non-trivial sufficient statistic for N. However 
an unbiased estimate is easily found, namely 


for 
k 
i=1 nj i 
] k 
=F > E[(N] = N. 


Using the approximate formula for the variance of N in the single 
sample case, (8), by a similar procedure to the derivation of (31) it is 
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found that 


In this inverse sampling procedure there are several parameters 
at the disposal of the experimenter —n, , k and s, , , 8. It 
may be desirable to choose them so as to minimize E(Li- » n;) while 
holding oy N fixed. This would achieve a fixed precision of estimation 
while minimizing the effort expended. However after arbitrarily fixing 
nu) and k, the optimum choice of the s; necessitates the solution of an 
algebraic equation of high degree (the degree increases rapidly with /:) 
which involves NV in a complex fashion. 

In lieu of general rules, Table III gives the properties of a number 
of simple designs. The approximate formula (32) was used to calculate 
ox/N. E(n;) was calculated recursively from the formula 


(NV + Vs; 
Em) — 8; + +1 


E(n;) = 


In order that this give a reasonable approximation it is necessary that 
s, be not too small. For this reason some of the more interesting cases 
with the initial s; very small are excluded. 

In the second model dealing with the inverse sampling procedure, 
at each stage n; — s; is predetermined. As before we write n; — s; = 6,. 
Thus 

(a) tin. — = 5; =0,1,--- ,k) 

(b) 8,82, M2, , m are random variables. 

Since ¢, is a parameter at each stage of the sampling procedure 
there is now independence between the successive random variables. 
Consequently the results determined in the single sample case are easily 
generalized to this model. Under the restrictions on n and 6 that (11) 
hold, it follows that 


+) 
(33) 
is an almost unbiased estimate. Also the variance is determined from 
(13) and the usual formula for the sum of the variances. Similarly 
k 
so) AD 

is derived from (15). 

It may be mentioned that a maximum likelihood estimate may be 
derived for this model which is analogous to that given by formula (19). 
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For again a telescoping occurs in the formula for the joint probabilities 
ping 
of the random variables viz. 


Prim, , Ms, *°* Ma) 


The maximum likelihood estimate of N is the solution of the equation 


where we write nm, for f, . In this case there is no simple sufficient 
statistic: the maximum likelihood estimate is a function of all the nx, . 
Since this is so and since the solution of (35) has a complicated distri- 
bution, the estimate -V seems more desirable. 


SEQUENTIAL TESTS FOR N 


The optimum sample census procedures evidently are sequential: 
furthermore the sequential procedure should permit a cheice of the 
design at any stage rather than merely a choice of whether or not to take 
further observations. For simplicity only standard sequential proce- 
dures are considered here, primarily in relation to tests for N. Such 
tests may be useful in control and management problems. 

For direct sample censuses with /(s,) small we consider 1; as pre- 
assigned and let / the number of stages in the census be a random vari- 
able. Since the s; have approximately a Poisson distribution, Wald’s 
theory [9] is applicable in a routine manner. 

In particular to test the hypothesis N against the alterna- 
tives VN > N, at a level of significance a, so that the power at a particular 
alternative NV, is L— 8 we proceed as follows: 


at any stage j after n, , , individuals have been sampled 
and , 8, , 8; tags recovered 
i 
(36) Reject it 5s; - 
log — log NV, 


(37) Accept if < Va__: 


7 log — log NV, 


Continue observations if $5), x, lies between the bounds in (36) 
and (37). 
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Here 


The bounds in (36) and (37) are themselves random variables: this 
is so because the s, are dependent in the model considered. However 
the conditional distribution of s; given s; (¢ < j) is of the Poisson form, 
approximately, with parameter n,t;/N. The setting up of the sequential 
probability ratio test does not require independence of the observations: 
it depends on the conditional distribution of the observed random 
variable at each stage. However the various optimum properties and 
associated results that Wald and others have proved for the sequential 
probability ratio test are not necessarily valid without independence. 
However approximate formulae could be determined for most situations 
from the formulae for the operating characteristic curve and for E(n) 
in the standard Poisson situation. Some of these may be found in [8]. 

In conclusion it may be reiterated that the sample censuses studied 
here have been limited to situations where the population is essentially 
stationary. This excludes populations where any substantial immigra- 
tion or emigration occurs. It does not necessarily exclude populations 
which are changing due to natural causes e.g. birth and death. If those 
born subsequently ‘to the initial tagging operation can be distinguished 
from the rest of the population, without undue effort and if the death 
rate is the same among the tagged and untagged groups, then the sample 
censuses conducted in one stage yield valid information on the popula- 
tion size at the time of first tagging. 

For suppose the probability of survival from time of tagging to 
time of sampling is p. Now consider the almost unbiased estimate in 
the direct single sample census model, where at the time of sampling (’, 
(N — t)’ = w’, N’ actually are surviving. 

Then 


(n + I(t + 1) t+1\_ 
= ar + 1 


-| t+1 


Since ¢’, u’ are independent and since it is reasonable to assume they 
have a binomial distribution 
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(38) 1+ ip (t+ 1) 


= N 


The denominator in the right hand side of (38) comes from the fact that, 
by direct summation 


for large t, and moderate p. . 

Similarly the unbiased estimate N = nt,s of the inverse sample 
census is unbiased whether or not mortality occurs in the population 
(provided tagged and untagged are proportionally affected). Even the 
variance of these estimates is only slightly modified unless the mortality 
is excessive. For example, the approximate formula (8) now becomes 


oi = N + 1) 
p s 


which is approximately N*/s unless p is very small. 

However if there is natural mortality the multiple sample census 
may be seriously affected—estimates of N in such a case may be mean- 
ingless. In fact this forms a basis for the possible estimation of natural 
mortality in an animal population—a fact utilized by Leslie and Chitty 
in a recent paper [10]. 

The references cited below are those referred to in the body or 
appendix of the paper. No attempt has been made ‘o compile a com- 
plete bibliography, even of the methodology in this field. In fact 
nowhere is such a bibliography available apparently; however the 
monographs of Ricker [14] and Schaeffer [16] give a large number of 
useful references. 
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APPENDIX 


I. To show the asymptotic normality of z as defined in formula 
(3) consider the moment generating function of z: 


By routine algebra and the usual manipulations it may be seen that 
In 1,(0) tends to 6/2 as s tends to infinity (e.g. ef. Cramer [5] pp. 
198-199). 

On the other hand the moment generating function of the random 
variable 2np is 


and In M,,,(0) = — s In [1 — 26 + R(p)] where > Oas p — 0. 
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lim = (1 — 20) 


which is the moment generating function of the x° distribution with 
2s degrees of freedom. 


Il. Maximum Likelihood Hquation. 
The maximum likelihood estimate is derived by setting the ratio 


where 


NIN —n, +58,) 


and also 


(42) 


so that S 1 according 
(43) 2 


also 


(41) q(N) S 1 


Thal 
Hye 
3 
f 
= 
S, R, 
: R, Rs 
{+ 
N 
S: 4 
N polis 
LS 
Rs 
cording as Vs, 


BIOMETRICS, DECEMBER 1952 


Let 


m = min ~~ M = max —— 


1 


In view of (40) and (44) the roots of (19) must lie between m and J. 
If one or more s; are zero JJ will not be finite. However unless all s, are 
zero (and hence S, = 0) it is evident from the form of equation (43) 
that g(N’) < 1 for sufficiently large N. S, = 0 may be neglected since 
Pr(S, = 0) will be extremely small. 

It is possible that. equation (19) has several real roots in the interval 
(m, M). Those for which q(N) is decreasing represent local maxima 
and one of these the extreme maximum. The large number of para- 
meters gives rise to a diverse number of possibilities but since in general 
N > >i... n; the following theorem will cover many cases that may arise. 
Theorem. If for N > m, R;/N**‘ is a monotone decreasing sequence 
and if 

k 
n; 


then the maximum likelihood estimate N* satisfies the following inequali- 
ties 


<2" 
Proof. 


since the left hand side is a weighted harmonie mean of terms 1,1, 5, 
and hence larger than the smallest term 


L= S, > 8, > m 
For 


so that V, > N* 
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TABLE 1 
NUMERICAL CoOAIPARISON GF CERTAIN INVERSE SAMPLE CENSUSES WITH s (NUM- 
BER OF TAGGED INDIVIDE BE RECAPTURED,) FIXED AND WITH 6 (NUMBER 
OF UNPAGGED INDIVIDUALS TO PE RECAPTURED) FIXED. 


2 3 > 
oN ox ot 
exact approximate 
100) 999} 101,000 3.25 BIZ! 
100, 4,950) 50 5.000 707.1 
| 500!) 475) 25 500; 5.0) woo!) .20 
| 100! 4,995! 5 | 5,000 2.2 226.1) 43 | .45 .55 
| 500 | 4,975 | 25 5,000 5.0 1000.0} 19 | .20 .20 
10° 1000 9,990! 10 10,000 | 3.2! 3162.3: 31.! | .35 
| 19,980 | 20 | 20,000 22 | 


1. 5 and s were chosen <9 that E(n) is the same (to the nearest integer) for both 
sampling plans, i.e., predetermining s or predetermining 6. 


; ot 

2. Calculated from the approximate formula ¢,, = VV 
Vs 
3. Calculated from the approximate formula ¢, = N - “7 


1. Calculated from formula (S) for the estimate where s is predetermined. 


(the approximation to formula (8)). 


6. Caleulated from formula (13) for the estimate where 4 ts predetermined. 


TABLE UU 
CUMULATIVE DISTRIBUTION OF X PoR VALUES OF 4, 5, 6. 7. 


P(N Sn 
0.0035 | 0.0590 0.3056 0.69044 0.94110 | 0.99685 | 1.0000 | 1.0000 
5 0.0012 0.0172 0.1052 | 0.3392 | 0.6608 | 0.8948 0982S 
6 - O.0001 0.0020 O.0627 | — .2010 
8.0001 8.0000 O825 0.1105 
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TABLE 
STANDARD ERROR OF THE ESTIMATE OF SAND TOPAL EFFORT REQUIRED 
BY VARIOUS MULTIPLE INVERSE SAMPLE CENSUSES. 


N = 10.000 N = 100,000 

Bs, Bs, | Ny E(n’)* oN E(n’)* oN 
| K(n') E(n') 

3 5,10,15,— —! 224 | 0.202 | 901 | 2.24 | 3385 | 0.60 
3 10, 10,10 — — 317 0.183 923 | 1.98 | 4020 0.46 
3 15,10, 5——j; 388 | 0.202 963 | 2.10 | 4595 | 0.44 
4 5,10,15,20—, 224 | 0.162 1131 1.43 | 3981 0.41 
4 12,12,13,13 —| 347 0.142 1155 1.23 { 4728 | 0.30 
4 20,15,10, 5 — 448 0.162 1208 1.34 | 5499 | 0.29 
5 5, 10, 15, 20,25. 224 0.135 1363 0.99 | 4617 | 0.29 
5 15,15, 15,15,15' 3880. 1394 | 0.83 5527 0.21 
5 25,20,15,10, 5; 500 | 0.135 | 1452 | 0.93 | 6363 | 0.21 
5 11,13, 15,17,19 332 0.118 | 1377 | 0.86 | 5170 | 0.23 
4 18,19,19,19 —, 425 | 0.116 1418 | 0.82 | 5794 | 0.20 
4 7, 15, 23,30 —' 265 0.134 | 1389 | 0.96 | 4835 | 0.28 
3 25,25,25 ——j| 500 0.115 1464 | 0.79 6373 
3 15, 25,35 — —j| 388 0.122 1425 | 0.86 5566 0.22 
2 37,383 —— —| 609 0.115 15389 0.75 | 7248 | 0.16 
2 25.50———! 500 ; 0.122 1513 | 0.81 | 6406 | 0.19 
1 7 —-—— —| 867 | 0.114 1732; 0.66 | 9508 | 0.12 


| 


* nn’ = total number of individuals sampled including those marked or tagged initially i.e. no. 

Since R; > 0 N, < Ni 

If N, < mthen NV, < N*. Consider the contrary case NV, > m and 
denote the other root of the quadratic equation (21) by N%. 


[the first inequality follows from the fact that 
27" <x 11 
Hence form < N < N, 


R, ite R 
Se < 0 und — (—]) > 4) 
Moreover 


provided only 2, < (2, 2)°. The right hand side of this equation is 
approximately of the order N/n, . 
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THE ESTIMATION OF RESPONSE-TIME DISTRIBUTIONS 
Il. MULTISTIMULUS DISTRIBUTIONS 


M. R. Samprorp 


Lectureship in. the Design and Analysis of Scientific Experiment, 
University of Oxford 


1. INTRODUCTION 


The first paper of this series (Sampford, 1952) consisted of an 
introductory discussion of the objects and methods of the analysis of 
response-time data. Among the topics held over for more detailed 
discussion was that of multi-stimulus distributions: these are dis- 
tributions occurring when a population or sample is subject to the 
action of two or more stimuli producing the same response, or producing 
different responses of which some are such that their occurrence pre- 
cludes that of the others. Such distributions are considered in this 
paper. 

2, GLOSSARY AND LIST OF SYMBOLS 


The previous paper contained a number of -definitions, some in- 
volving symbols and some particular uses of words, which will be 
followed consistently throughout the series. For convenience, these 
definitions are repeated here, together with others occurring for the 
first time elsewhere in this paper. The symbols listed here are those 
which are used consistently in more than one context, and which are 
neither completely standard nor self-explanatory. 

Glossary. 
The following words are used in a specialised sense: 


Accident: An animal .. said to be accidentally killed if its death is 
caused by the experimenter or by the conditions of the experiment, 
but not by the treatment, or if it is injured in the course of the 
experiment and is destroyed. Thus deaths from such causes as 
post-operative shock and the after-effects of anaesthesia would be 
classified as ‘accidents’, as well as obvious accidents in, say, in- 
jection or feeding. 

Initial event, Reaction, Response-time: The initial erent is the origin 
from which time is measured. It is usually the time of treatment, 
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but it may be some appropriate later time such as the time of 
appearance of a symptom other than the reaction. The reaction is 
the intended effect of the treatment, and the response-time is the 
time elapsing between the initial event and the reaction (Sampford, 
1952: § 2.1). 

Mortality, Kill, Survivors, Survival-time: Many response-time problems 
and methods, even those of general applicability, are most con- 
veniently described by assuming the reaction to be death, and by 
considering the time-mortality distribution, for which the response- 
time is the survival-time, while the percentage kill and percentage 
surviving at any time may be interpreted as the percentage which 
have or have not responded at that time. 

Sacrifice: An animal is said to be sacrificed if, for any reason, it is 
deliberately removed from the experiment before showing the 
reaction. 

(Absolute) survival, Immunity, Recovery: UHere again, although the 
terminology of time-mortality is used, these concepts are equally 
applicable to other responses. As indicated above, all individuals 
which have not shown the reaction at a particular time are said to 
survive to that time; for certain treatments, however, there may be 
some individuals that fail to respond, even though observation is 
continued so long after the last recorded reaction that the prob- 
ability of any further response is negligible. Such individuals are 
said to survive absolutely; those unaffected by the treatment are 
said to be immune, those affected, but not enough to show the 
reaction, are said to recover. 

Truncation, Partial Truncation: A series of data is said to be truncated 
if individual observations are available only over part of the range, 
so that the information available on any individual can take one of 
the three forms: 


T 
r= 7, 
r> 


where 7', and 7’, are the same for all individuals in the sample. In 
particular cases T, may be the lower, or 7 the upper, limit of 1, 
so that the data are truncated at one extreme only. In time-response 
work the numbers of observations in the tails (i.e. for which 7 < Ty , 
or 7 > 7;,) are usually known, so that the data are of the type which 
Hald (1952a; § 6.9) has called censored. When the points of trunca- 
tion vary from individual to individual, the data are said to be 
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partially truncated. Wor example, the only information about. the 
survival-times of two individuals A and B might be that they were 
greater than r4 and rp, respectively, where 74 need not be equal to rz . 


List of standard symbols 


7 ~~ an observed response-time, us defined above. 
—- a time-metameter, assumed normally distributed in 
the absence of interference. Usually ¢ = 7, log 7, or 
1/7. It is assumed, where relevant, that ¢ is an in- 
creasing function of 7, and some of the methods here 
derived depend, in detail, on that assumption. When 
{is a decreasing function of 7, such as the reciprocal, 
the sample t-values should be subtracted from zero or 
another suitably chosen constant and analysed by the 
methods as given. 
u, o — mean and standard deviation of the distribution of a 
time-metameter, usually but not necessarily normal. 
» — standardized normal deviate, defined by 


i - 
=a+ Bt 
o 


a, 8 — alternative parameters for normal distribution, de- 
fined as above, so that 


a = —p/o, B = I/e 


m, s — estimates of py, 
a, b —- estimates of a, 8 
y(t) — estimate of n(t): 


m 


y(t) = ‘= or a+ bt. 


Z(n) — ordinate of standard normal distribution correspond- 
ing to deviation 7: 


Un) = 


V 


P(m), Q(n) --- left and right tails of normal distribution, measured 
from 7: 
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(When this relation is used to define 7 as a function of P, q is ealled the 
normal equivalent deviate, or N..D., of P.) 


o. AD 
v(m) dy log Q Q(n) 

Mn) = log Q = v(m) — 


= v(n) — 
The func’ions », A, and ¢ are tabulated in the appendix to this paper 
(—5.0(0.1) —3.0(0.01)3.0(0.1)5.0). Hald (1952b; Table X) has tabu- 
lated a function 
= (—3.0(0.1)2.0). 


Cornfield and Mantel (1950) have tabulated related functions for use 
in probit analysis. They tabulate (against probits; 1.0(0.01)9.0) the 
functions A(x) and — w(x) for maximum and minimum response, 
where 


A(x) 
Minimum 
—w 


v(x — 5) 


3. THE SIMULTANEOUS OPERATION OF TWO STIMULI 


— 5) 


Only distributions resulting from the action of two stimuli will be 
considered; the extension to greater numbers is usually obvious. The 
word ‘stimulus’ is here used in its widest sense, as including not only 
applied treatments, but the complex interactions of circumstances 
which give rise to accidental or natuial deaths. 

There is an important distinction between distributions in which 
the particular stimulus responsible for the reaction of any given in- 
dividual is known, and those in which it cannot be determined. The 
first group includes, of course, all those in which the stimuli give rise 
to different reactions: for example, an animal treated with carcinogen 
may develop a tumour, or it may die from natural causes before a tumour 
appears. This group alsv includes those cases when either stimulus 
can kill the subject, but the cause of death is apparent: in particular, 
that in which the two causes of death are an applied treatment and 
accident. The second type of distribution, for which the cause of death 
cannot be identified, occurs, for example, when a group of animals is 
treated with a mixture of two materials, either of which is capable of 
inducing convulsions, or when an experiment on the effects of atomic 
radiation continues for a period so long that the experimenter is unable 
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to say whether any particular animal dies as a result of the treatment 
ov from natural causes. ‘The problems of estimation appropriate to 
distributions of the first type are usually less difficult, and give rise to 
considerably Jess laborious calculations than those connected with the 
second type. The two types are exemplified in this paper by the 
‘accidental death’ model of § 5 and the ‘natural death’ model of § 6. 

In neither case does the form of the analysis depend at all on the 
nature of the reactions concerned. In view of this, and for convenience 
of terminology, all methods developed in this paper will be discussed 
in terms of the typical reaction of death. 


4. THE JOINT OPERATION OF TWO CAUSES OF DEATH 
4.1 Types of joint action 


When two causes of death operate simultaneously on a population 
the resultant distribution of survival-times depends on the type of 
joint action in operation. The possible types of joint action have 
been discussed by Bliss (1939) and others, in connection with quantal 
responses to mixtures of poisons. I shall follow the nomenclature of 
Plackett and Hewlett (1952), who distinguish four types, classifying 
the action as similar or dissimilar according as the two poisons affect 
the same or different physiological systems, and as simple or complex 
according as the effects of the poisons are or are not additive. Their 
four types are thus: 


(i) Simple similar action: the two poisons have a common site of 
action; if, too, the application of a dose 2, of poison 1 produces 
an effective dose u, of the active principle at the site of action, 
and a dose x. of poison 2 produces an c¢ffective dose u, , the 
simultaneous application of doses 2, and a, of poisons 1 and 2 
produces an effective dose u, + us. . 

(ii) Complex similar action: there is again a common site of action, 
but the effective dose at the site of action is no longer u, + wu, ; 
it may be greater (synergism) or less (antagonism). 

(iii) Simple dissimilar (or independent) action: the two poisons 
operate at different sites; furthermore, the presence of one 
poison affects neither the amount of the other reaching its site 
of action, nor the magnitude of the response at the site. 

(iv) Complex dissimilar (or dependent) action: there are again 
different sites of action, but the two poisons interact. One or 
both may inhibit or intensify the action of the other, or poison 
1 may inhibit poison 2, while poison 2 intensifies the effect of 
poison 1. 
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It remains to consider the types of distribution of survival-times to be 
expected from these four types of joint action. 


4.2 Similar action 


In the field of response-time investigation, similar action seems 
likely to occur only when the two treatments are doses of two poisons 
having a similar or identical active constituent. In this case the dis- 
tribution of response-times would be simply that corresponding to the 
effective combined dose; thus the three response-time distributions 
resulting from the separate and simultaneous application of the two 
poisons would correspond simply to three different effective doses of the 
same active principle, and as such might be expected to have similar 
forms. In particular, if the two poisons applied individually give 
distributions normal in some time-metameter ¢, the distribution result- 
ing from their joint action would probably also be normal in the same 
t, the three distributions thus differing only in their means and, possibly, 
their variances. It may be assumed that the experimenter is always 
free to apply the two poisons separately and together at will; this situa- 
tion thus presents no peculiar problems of estimation. The response- 
time may be used, in the same way as any other quantitative response, 
to test for interaction between the poisons, but a discussion of such tests 
is more conveniently postponed until after a consideration of the re- 
lationship between dose and response-time. The remainder of this 
paper, therefore, will be confined to the discussion of problems of 
dissimilar action. 


4.3 Dissimilar action. 


In a population subject to two causes of death having different 
sites of action, each individual may be considered as having ‘potential’ 
survival-times ¢, and ¢, corresponding to the two causes, where ¢, and 
t, are measured from a common origin, and are jointly distributed in 
the population according to the frequency function 


S(t, dt, dt, (1) 
If cause 7 operates alone, survival times are distributed in the form 
dt; (2) 


where the f, are marginal distributions of (1) if and only if the joint 
action is independent (i.e. of type (iii) as defined in § 4.1). When both 
causes operate simultaneously, only one of ¢, and é can be realised, 
and death results from cause 1 or cause 2 according as t, < ¢, ort, > l 
respectively. Thus an individual dies from cause | in t — (t + dt) if 
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{, lies in that range, and provided it has not atready died from cause 2: 
that is, provided 4, > ¢. The probability of this situation is 


Similarly, the probability of a death from cause 2 in the same period is 


an} ae (3b) 


and the distribution function of survival-times for deaths from both 
causes is 


su, dt at (3c) 


In the particular ease when ¢, and f, are independent and uncorrelated 
ty) = S(t) fo(te) (4) 
and the functions (3a), (3b), and (8c) become 


— fio) de 
[1 — de (5) 


and H(t) = [Fi + FAO — 


where /’;(t) is the distribution function corresponding to the frequency 
function 

If the distribution functions F,(¢;) were known, it would be theo- 
retically possible to devise appropriate tests to distinguish between the 
various types of joint action. This might be of some practical use 
when the causes of death are two poisons, which can be applied ex- 
perimentally, first singly, and then together. A possible test for this 
type of experiment is suggested in § 6.11, for the case in which the dis- 
tributions F,(¢;) are normal. More usually, however, the separate 
distributions are unknown, and the object of the analysis is to estimate 
one of them, which we may conveniently assume to be F,(t), from 
data in which its operation has been disturbed by the simultaneous 
operation of a second cause of death, having the ‘nuisance’ distribution 
of survival-times F’,(¢,). When the bivariate p.d.f. f(t; , &) can be 
satisfactorily defined in terms of the parameters of {,(é,), these param- 
eters can theoretically be estimated from the sample; this is so, for 
example, when the joint action of the stimuli is independent and un- 
correlated, so that f(t,, 42) has the form (4), or in the particular ease of 
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independent action in which 4) is bivariate normal, defined in 
terms of the means and standard deviations of (4) and f.(é:) and a 
further unknown parameter p. The adequacy of the estimator in any 
particular case will, of course, depend on the size of the sample and the 
number of other parameters involved. In the case of dependent action, 
however, an attempt to express /(f, , /:) in terms of the parameters of 
fi(4) would probably involve various complicated assumptions, very 
difficult to justify, on the nature of the joint action of the stimuli. 
This limitation is not, in practice, as serious as it appears in theory. 
The three most common ‘nuisance’ causes of death are disease, accident 
and sacrifice, and natural causes. In the case of disease, admittedly, 
the difficulties of analysis are considerable. The effects of the disease 
will almost certainly hasten the effects of the treatment, and an attempt 
to estimate /’,(¢,) is unlikely to sueceed. If the disease is epidemic and 
makes its appearance late in the experiment, the data may be truncated 
at the time of onset, and F',(¢,;) estimated from the truncated data alone; 
in other cases the only course open to the experimenter seems to be a 
repetition of the experiment using healthy stock. (This assumes the 
existence of a healthy stock: if the disease is ineradicable it must be 
regarded as a natural cause of death.) 

Deaths as a result of accident during treatment or inspection may 
occur in any experiment which lasts longer than a few days, and animals 
are sometimes sacrificed while an experiment is in progress, either to 
reduce the size of the experiment or for autopsy. In almost all such 
cases the joint action will be independent, and provided reasonable 
precautions have been observed by the experimenter survival-times for 
the two causes should be uncorrelated. (Correlation might oecur, for 
example, as a result of cannibalism, which is fairly common among 
experimental insects. The victims would tend to be the weakest 
animals, or those most affected by the treatment, and these might be 
expected to have the shortest potential survival-times. Neglect of this 
fact would result in an overestimation of the mean treatment survival- 
Bae time. Similar remarks apply to deaths occurring as a result of fights 
4 when experimental animals are kept more than one to a cage). Animals 
p to be sacrificed should, of course, be chosen at random from all survivors 
at the time of sacrifice. 

In long-term experiments, animals may die from natural causes. 

Here again, the joint action is likely to be dependent, with the operation 
: of the treatment accelerating natural death, or, at best, independent 
correlated, as the weakest animals, which would die first if untreated, 
will probably succumb most readily to the treatment. [stimation of 
F,(t,) will then be difficult in the latter case and practically impossible 
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in the former. This difficulty, however, is not as serious as if may 
appear at first. The distribution /\(4,) refers to survival-times in the 
presence of cause | only, cause 2 being absent. When ‘cause’ 2 is the 
natural mortality among the population, /,(¢,), so defined, can never 
be realised in practice. The treatment cannot but operate in a natural, 
and hence mortal, population; what is required is not the hypothetical 
’,(t,), but the distribution of deaths superimposed by the treatment 
on the natural mortality distribution of the population. In practice, 
the actual cause of death is often uncertain, so that the data consist of a 
single sample distribution estimating f(t, , 4). 

It therefore seems reasonable to attempt to approximate to f(d, , 6) 
by a bivariate distribution of the form 

Wty) fo(te) 

where f.(t) is the p.d.f. of survival-times in an untreated population 
and g(é,) the ‘superimposed’ distribution function, and to attempt to 
estimate g(¢,) rather than the hypothetical /,(é,). The suecess of this 
approach would depend on the existence of a fairly simple form of 
g(t;) providing an adequate approximation. This problem is discussed 
in more detail in § 6, with reference to the case when /,(é,) is normal 
and /(t; , t2) is bivariate normal; in this case an assumption of normality 
for g(4,) appears to give a mathematically adequate model. 

The two most important situations, then, are those which may be 
described as the ‘accidental death’ and the ‘natural death’ situations. 
These will be discussed in the two following sections. 


5. THE ‘ACCIDENTAL DEATH’ MODEL 
o.L Statement of the problem. 


The situation discussed in this section is as follows: a distribution of 
survival-times to a particular treatment, assumed normal in terms of an 
appropriate time-metameter, is modified by the occurrence of death 
from another cause. The two causes of death are distinguishable, and 
the joint action is assumed to be independent and the survival-times 
uncorrelated. The most common ‘nuisance’ causes giving rise to this 
situation are accident and sacrifice. 


0.2 Estimation of parameters. 


It is assumed that the distribution of treatment survival-times is 


h(t) dt = — = BZ(y) dt 


= — Qi). 
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The distribution of survival-times from accident or some other 
secondary cause is 


f2(t) dt 


and is assumed not to depend on the paramters a and 8. Then from 
formulae (4), since the joint action is assumed independent and un- 
correlated, 


likelihood of a death from treatment at time ¢ = [1 — F,.(t)]8Z(n) 
likelihood of a death from secondary cause at time t = f,(t)Q(n). 


Then if N individuals are treated initially, and, of these, n die as a 
result of the treatment at times ¢; , and n’ = N — n die as a result of 
the secondary cause (e.g. are accidentally killed or sacrificed) at times 
t; , the likelihood of the sample is 


L 


— {TT £20000) 


TT 2m) A 
where A(t) is independent of a and 8, so that 


log L = n log 8 — + log Qa) + BO 


Now = a+ £t; dn/da = 1; = t. 
The functions »(), A(n), and {(m) are defined as follows: 


v(n) = Z(n)/Q(n) = —dlog Q/dy 
= ¥°(m) — = dv/dn (6) 


= v(m) — nA(n) 


They are of importance in this and subsequent sections, and are 
tabulated in an appendix. References to tables of related functions are 
given in § 2; some properties of v(m) and \(n) have been discussed else- 
where (Sampford, 1952). 

The maximum likelihood equations for a and 8 are 


t’ 
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and are independent of the distribution F,(). They cannot be solved 
explicitly, and it is necessary to use an iterative method. 
The Maclaurin expansion gives, to the first order of small quantities: 


log L . log L log 
(8) 
da 0B 


where the coefficients of the increments and the right-hand members 
may be evaluated in terms of a first approximation a, , b,; and used to 


calculate increments leading to a, , b, . 
Now 
poe 
log L 


| 
| 


where \ is the function defined above, so that equations (7) take the 
form 


{n + ae 6a, + 6b, 
a >» 


t hd t q (9) 


b, J 


where y, = a, + is an estimate of and are functions of , 
and da, , 6b, are the increments to be added to a, , b, to yield the revised 
estimates a, , . 

Introducing a general variable ¢ to inelude both ¢ and ¢’, a weighting 
function w, such that 

w = 1 for an individual dying from the effects of treatment 
= \, for an individual dving from the secondary cause, 

and a general summation >> operating over all ¢ and 1, equations (9) 
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reduce to 


6a, w+ db, v= wy, + > AW vy 
ba, >> wt + db, we} (10) 


Adding to both sides of the first equation )> wy, , and to the second 
equation n/b, + >> wy:t, it follows that 


a, > wt = — > — dy} 
(11) 
and finally, writing 
dwt 
=n Am = Cy) 
and 
y =a’+ W(t — 2) 
so that a’ = a + bt 
equations (11) reduce to 
a; = 
(12) 
2n 
= Se — + 


If necessary, further cycles can be carried out, until stable values 
of a’ and b are obtained. The asymptotic variances of a’ and b are 
then estimated by the reciprocals of the denominators in (12), as 
calculated for the final cycle. 

For practical calculation, it is convenient to return to a stage rather 
earlier than the final form. 

tand are the same for all cycles, and are calculated once 
for all at the beginning of the analysis. 

For each cycle it is necessary to tabulate, for the animals acci- 
dentally killed, y, A, and ¢. From these values, >>" A, 35" At’, and 


ae 
S 
f 


RESPONSE-TIME DISTRIBUTIONS 319 


are caleulated, and added ton, to give wt. 
wl. and et’ are also calculated from the table. Finally 
t, wt — 1°, — 0, 2n/b, and n/b; are calculated, and equa- 
tions (12) used to calculate a3 and b. . 


5.3. Estimation of secondary distributions. 


The method of § 5.2 assumes only a general form for the distribution 
f(t), and is thus general in its application, provided only that. the 
joint action of the two causes of death is independent and uncorrelated. 
It is not even necessary that the form /,(é) should be identical for all 
individuals. There are a number of important special cases. 

(i) Truncation. Here no observations are made after time dy , so 
that all surviving animals may be considered as being sacrificed at 
that time. Here 


= 1; =0, 


Thus the method of § 5.2 provides an alternative to the various com- 
puting routines already developed by other writers (Bliss & Stevens, 
1937; Hald, 1952a; Cohen, 1950) for the caleulation of maximum 
likelihood estimates of the parameters of a normal distribution truncated 
at an upper limit. For data truncated at the lower end of the scale 
some slight modification is necessary. This problem will be discussed 
in more detail, and the relative merits of the various methods will be 
considered, in the third paper of this series. 

(ii) ‘Partial truncation’. This situation, corresponding to the 
‘incomplete graded response’ problem described by Ipsen (1949), 
oceurs When the distribution functions of different sample members are 
truncated at different points. Thus the only information available 
about two members A and B might be that they were still alive at 
times 7, and 7, , where 7, is not necessarily equal to rt, . In this ease 
J,(t) may be defined for each individual 7, as 


f(t) = 1; = 0, 


and the method of § 5.2 may be used to estimate the treatment mortality 
distribution, provided ¢; is uncorrelated with the expected treatment 
survival-time in the population. 

Partial truncation frequently occurs when response-times are 
measured not from the time of treatment, but from some event which 
may occur at any time after treatment. If the experiment is terminated 
artificially, individuals that have still to show the reaction will have 
histories starting from different initial events, so that when all survival- 
times are reduced to a common origin, the data will appear partially 
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truncated. For example, if animals are treated with a carcinogen, 
tumours will appear at various times after treatment. Hf the ‘initial 
event’ is the appearance of a tumour, and the respouse-time is its 
survival-time to disappearance or ‘regression’, any tumours still sur- 
viving when observation is suspended will date from different times of 
appearance, so that the sample of survival-times will be partially trun- 
cated. The methods of this section can, of course, be used in such a 
case only if the experimenter has good reason to believe that survival- 
time is in no way dependent upon date of appearance. | 
(ili) Sacrifice. Several different situations are possible. The experi- 
menter may decide in advance to sacrifice animal 7 at time 4; , where 
the {; may be all the same (truncation) or different (partial truncation), 
and where some of the ¢; may be infinite, so that observation must be 
continued until these animals have died from the treatment. Alter- 
natively he may choose to sacrifice proportions p, , or numbers n,; , of the 
survivors at one or more times {, . Provided all choices are made at 
random, the hypothesis of uncorrelated independent action will hold. 
Formal expressions for f.(¢) can be written down for all these models, 
but there seems no reason to do so, as in no case is any question of 
estimation involved. 
‘iv) Accident. If a knowledge of the accident distribution is re- 
quired, and if the experimenter is prepared to hazard a guess at the 
form of the distribution, depending on one or more unknown parameters, 
the functions f,(t), F.(t) may be replaced, in the development of § 5.2, 
by their supposed forms. The functions A(t), B(é) will then depend on 
the parameters of the accident distribution, and may be maximised to 
yield estimates of these parameters. 
The most plausible assumption of a form for f.(¢) seems to be that 
any individual which survives to time 7 has a probability «x dr of being 
accidentally killed in the interval, 7 — 7 + dr. This gives 


Sr) dr = dr 
0< 
1- Fr) =e" 


= Fr) =0 7 <0 
Inserting these values in the likelihood function of § 5.2 
A(t) = 


where >> 7 is a summation over ¢ and (’ of their values in true time, if 
= ‘log r 
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Maximizing A(é) for x gives the maximum likelihood estimate 
anf 
(13) 
T 
Also 
L _ log A _ 
On 
so that the asymptotic variance of & is 
= (14) 


(v) Second cause of death normally distributed in t. In this case, 
provided the joint action is independent and uncorrelated, the method 
of § 5.2 may be used to estimate F',(t). If the form of F.(#) is of any 
interest, it too may be defined in terms of parameters a, and 8, , which 
then appear in the function A(t). The method of § 5.2 may then be 
applied to A(t) to obtain maximum likelihood estimates of a, and B, . 

(vi) More than two causes of death present. If three or more causes 
of death are present, provided that all secondary causes are distinguish- 
able from the main cause, and that their joint actions with it are in- 
dependent and uncorrelated, they may be considered as a single cause 
for the purpose of estimating F,(é), and the distribution F.(é) will then 
be the distribution produced by the joint action of all secondary causes. 
If the secondary causes are themselves independent, uncorrelated, and 
distinguishable, the results of § 4.2 may be extended. For example, if 
there are three causes of death, the formulae (5) become: 


[1 — FACIL — at 


[1 — — at 
(15) 
(1 — — dt 


and 


= 1 — [1 — — FAO) — 


where the first three results are the probabilities of death in the interval 
t — t + dt from causes 1, 2, and 3 respectively, and H(t) is the joint 
distribution function of the three causes acting together. The method 
of § 5.2 may then be extended to estimate F,(é), including in the n’ 
animals dying at times ¢’ all animals dying from all causes other than 
cause 1, The distributions /’.(f), F;(€) may also be estimated if required. 
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For example, suppose n animals die of cause 1 at times ¢; , n’ as a result 
of accident (cause 2) at times ¢’ , and n” are sacrificed at times t/’.. Then 
if cause 1 has the normal distribution of § 5.2, and cause 2 the accident 
distribution of § 5.3 (iv), the equations (12) of iteration for a and 8 are 
unchanged, except that the times ?’’ are also included in any summation 
involving ¢’, and the equations (13) and (14) for % and V(x) are also 
unchanged, except that ) 2 7 now represents a summation over the 
true times corresponding to /, ¢’, and ¢’”’; if t = log r. 


De + De’ + 
5.4. Tests of goodness of fit. : 


If the estimation procedure has included the estimation of the 
secondary distribution F(t), the joint distribution function H(é) may 
be built up from formulae (5). In general it will not be possible to 
calculate exactly the expected numbers dying from the separate causes 
in any finite interval, but the expected number dying from both causes 
together can easily be calculated, and compared with the observed 
values by a x’ test in the usual way, degrees of freedom being obtained 
by subtracting from the number of intervals one more than the total 
number of parameters estimated in F,(f) and F,(t). This method, how- 
ever, provides only a rather insensitive test of goodness of fit of F,(t) 
if « fairly large proportion of total deaths are due to secondary causes, 
and is inapplicable when some of these are due to sacrifice. An alterna- 
tive test, however, may be derived as follows. The probability x; that 
an animal surviving at the beginning of any interval will die during the 
interval, ignoring the possibility of deaths from secondary causes, can 
be calculated from the estimated F,(i). If, then, the time-scale is 
divided into intervals, starting a new interval whenever a death from 
secondary causes occurs, the number n,; surviving at the beginning of 
each interval is known, so that the expected number of deaths occurring 
from cause | in each interval can be calculated as n;x, , and compared 
with the observed number, /;. In practice, unless the secondary deaths 
are infrequent, this is likely to lead to very small expected numbers 
n,x,; , so that it is more convenient to choose the intervals of adequate 
length, calculating the expected numbers dying from cause | as n/n, , 
where n; , the average number at risk in the interval, is given by 


ni = n; — one half the number of individuals dying from 
secondary causes in the interval. (16) 


This is, of course, an approximation, but one which is not likely to be 
serious, provided care is taken to choose any times at which considerable 


al 


a 


4 
f 
«ORI 
j 
sa 
be 
ou 
in 
i 
= 
4 
in 
{ 
i 
: 
‘ 
q 
are 
| 
| 
| 
a 


— 


RESPONSE-TIME DISTRIBUTIONS 323 


sacrifices are made as end-points of intervals. (In fact, intervals should 
be chosen so that n’x,; (1 — x;) isnot toosmall.) The test is then carried 
out as follows. The whole range of the time-scale is divided into k (say) 
intervals. For each interval, the a priori probability of a death from 
‘ause 1, , is calculated x; = 1). Then 


and n{ is calculated as above. -Then the number /, dying from cause 1 
in the 7th interval is binomially distributed, with approximate expected 
value n,x; , so that the quantity 


Xx = 


nin — x;) (I) 


is approximately distributed as x’, with degrees of freedom (k — m — 1), 
where m is the number of parameters estimated from the sample which 
have been used in calculating the z; . In particular, if the distribution 
F(t) is assumed normal, and parameters a and 8 have been estimated 
from the sample, the number of degrees of freedom is 3 fewer than the 
number of intervals. Very small values of n/x; (1 — x;) are undesirable; 
intervals should be chosen to avoid such values, and large contributions 
to (17) from intervals with nix;(1 — 2,) less than, say, 5, should be 
accepted with caution. (It will be noticed that.the Ath interval, with 
upper limit 7 = ©, provides no contribution to (17); all animals sur- 
viving at the beginning of the interval, and not killed by secondary 
causes, must die from the treatment during the interval, so that both 
the observed and the expected number killed are equal to 1, . Alter- 
natively, all surviving animals may have been killed at the end of the 
(k — 1)th interval, so that both nj and l, are zero. In either case the 
term of the summand in (17) is identically zero; the last interval should, 
however, always be included when counting the intervals, to give the 
correct number of degrees of freedom.) 


5.5. Example. 


The data of table I consist of the number of shrimps (Lysmata 
seticaudata) moulting, and dying before moulting, on successive days 
after treatment with a hormone preparation. Observation was con- 
tinued until the end of the ninth day, at which time all survivors may 
be taken to have died. Consideration of other similar series suggests 
that after such treatment the response-times are log-normally dist rib- 
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uted. All animals dying on any one day are taken as dying at the 
midpoint of the interval, and the time metameter is taken as the log- 
arithm of the midpoint, as recommended in the first paper of this 
series, rather than as the mean of the logarithms of the endpoints. All 
values of t have been increased by 1 in order to avoid negative values. 
(As these calculations are included simply as an example of the methods 
of this section, the recommendation of the previous paper, that all 
animals responding in the interval 0-1 should be considered as falling 
in a truncated lower tail on the log-scale, has been disregarded, and 
such animals are assumed to have an R.T.M. of (1 + log 0.5). Further, 
as there is considerable doubt whether the required assumption of 
independent and uncorrelated action between cause of death and cause 
of moult is justified in practice, these calculations should be regarded 
as an example of method only.) When, as here, the numbers accidentally 


TABLE I 


NUMBER OF SHRIMPS MOULTING AND DYING WITHOUT MOULTING ON SUCCES- 
SIVE DAYS AFTER TREATMENT WITH HORMONE PREPARATION. 


a Day No. moulting No. dying T t 

1 8 5 0.5 0.70 
2 7 4 1.5 1.18 
3 4 1 2.5 1.40 
4 4 6 3.5 1.54 
5 8 3 4.5 1.65 
6 7 9 5.5 1.74 
7 5 9 6.5 1.81 
8 2 6 7.5 1.88 
9 2 5 8.5 1.93 

Survivors 34 9.0 1.95 


killed are large, provisional approximations may be obtained as shown 
in the first part of Table III. For each interval the average number at 
risk, n{ , is calculated as the number surviving to the beginning of the 
interval, less one-half the number killed in the interval. From these 
values and the observed numbers moulting 1; , are calculated 
(i) the proportion of ni moulting in the interval, p{ 
(ii) this quantity expressed as a proportion of the original popula- 
tion, 
Pi = and 

(iii) the proportion surviving at the end of the interval, assuming 

_ zero death-rate 


ae 
it 
| 
| 


it 
1e 


RESPONSH-TIME DISTRIBUTIONS 


325 


(p; , p; , and Q/ are the empirical values of the quantities which appear 
in the second half of table Il] as x; , x, , and Q, , there calculated from 


the estimated parameters.) 


The N.E.D’s y; corresponding to the Q/ can then be plotted against 
the values of ¢ corresponding to the ends of the intervals, and an empirical 


line drawn by eye to provide provisional estimates of a and 8. In this 
case the provisional values are 
a, = —3.21, b, = 1.57 
The calculation then proceeds as in Table II. 
TABLE II 
ANALYSIS OF DATA FROM TABLE 1. 
FOR INDIVIDUALS MOULTING: 
n = 47 Xt = 67.67 De = 104.8655 
FOR INDIVIDUALS DYING: 
t’ ni y nir nini’ nit 
0.70 5 —2.87 0.472. 0.330 1.215 
1.18 4 —1.96 1.063 1.254 2.154 
1.40 1 —1.54 0.367 0.514 0.655 
1.54 6 —1.27 2.592 3.992 4.279 
1.65 3 —1.06 1.442 2.380 2.243 
1.74 9 —0.89 4.673 8.131 6.918 
1.81 9 —0.76 4.934 8.931 7.030 
1.88 6 —0.63 3.457 » 6.498 4.739 
1.93 5 —0.53 2.978 5.747 3.971 
1.95 34 —0.50 20.489 39.954 27.040 
First approximation a,= —3.21 b= 1.57 
42.467 > w= 89.467 nig= 60.244 
> vat’ = 77.731 > wt= 145.401 nigt= 108.6628 


> nidt’? = 144.0363 > wt? =248.9018 


t= 1.625 w(t—t?= 12.5973 — n’¢(t—t) = — 10.7663 


n/bi= 19.0677 2n/b,= 
Total 31.6650 
b= 1.55 V(b.) = 1/31.665 = 0.031581 
a, = —0.673 = V(a2) = 1/89.467 = 0.011177 
a, = —3.19 V(a,) = V(at) + f(b) = 0.094571 


59.8726 


49.1063 
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As the values of a and 6 are practically unchanged by the first cycle, 
no further cycles are necessary. The test of goodness of fit is shown in 
the second part of Table III. The last two days are taken as a single 
interval to avoid very small values of n;x; (1 — 2,), the number at 
risk being taken as (number surviving to beginning of eighth day, i.e. 
to 7 = 7) — 4 (number dying in eighth and ninth days). The x’ 
of 6.21, with 6 degrees of freedom, suggests a degree of goodness of fit 
which is to some extent belied by the systematic nature of the deviations 
from expectation. (However, comparison with other similar series 
suggests that this is a purely chance effect; over all series the assump- 
tion of normality for log + provided a very satisfactory fit.) 


6. THE ‘NATURAL DEATH’ MODEL — TWO NORMAL DISTRIBUTIONS 
6.1. Introduction. 


In this section a distribution of survival-times after some treatment 
is assumed to be modified by the presence of another cause of death, 
such that the joint distribution of ‘potential survival-times’ of § 4.2 is 
bivariate normal in terms of some time-metameter. Where the two 
causes of death are distinguishable, the data may be analysed by the 
method of § 5.3(v) if the survival-times ¢, and ¢, are uncorrelated; other- 
wise by an extension of the methods of section 5. If the causes of death 
are indistinguishable, however, the data consist only of a sample from 
the multi-stimulus distribution resulting from the joint action; a dis- 
tribution depending on five parameters (or four, if the survival-times 
ure assumed uncorrelated). To estimate even four parameters from a 
single sample of practicable size would be difficult, but the situation is 
eased by the fact that in the most important case, that in which the 
secondary deaths occur as a result of natural causes, a control series will 
usually have been run alongside the treated, and may be used to provide 
estimates of the parameters of the marginal distribution /,(t,). Before 
discussing the estimation procedure some knowledge of the form of the 
multi-stimulus distribution, and particularly of its moments, is desirable. 
These are discussed in the following section. 


6.2. The moments of the resultant distribution. 
In the notation of § 4.2. 


and the resultant distribution of all deaths, 
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h(t) dt = {h,(t) + dt 
where 


=> 


— p 


and h,(t) is the corresponding function obtained by interchanging the 
suffices 1 and 2. 
The moment generating function of H(t) is given by 


= + ¢2(6) 


where 
= dt. 
Then 
1 
(0) = af at 
2(1 — p) 


which, writing v = u — ¢ and reversing the order of integration, reduces 
after considerable manipulation to 


292 4 6 
g)- = 6 Me My (o; 
Voi — + 


and similarly 


Voi — + 
Thus 


Voi — 2p0,0, +- 
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where 


Be — 


Differentiating with respect to 6, and setting @ = 0, we obtain the first 
four moments of //(t): 


= miP(A) + — Voi — + 03 Z(A) (193) 
uss = + mi)P(A) + (02 + Q(A) 

— (ur + m2) Voi — + 03 (193i) 
= (Buoy + + + 


— (ui + + Vor — + ZA) 


2oi — + (1 — — + 203 
Voi — + 03 


uss = (301 + + wi)P(A) + (302 + + u2)Q(A) 
— (ui + + + Voi — 2po.02 + Z(A) 
_ — + 3o2)u, + (30; — + A 
Vo; = + ( ) 


oF p ){(or — + Cz P7102) Mo} 
— + 02)” Z(A) (19iv) 


(19iii) 


6.3. Estimation by moments. 


A complete estimation of all the parameters of the distribution 
H(t) from the sample moments of the joint distribution is impracticable. 
If, however, we know, or are prepared to guess, the value of p, and 
#2 and o2 can be adequately estimated from a control sample, the first 
two moments may be used in a simple iterative method for estimating 
nu, and o,. For from (19i) above 


Mis — + Vo; - - +o o> Z(A) 
Mi = + P(A) (20) 


and from (19i) and (19ii) 
(oi + wi)P(A) + (63 + 

— + me) + — 
+ o:Q(A) + + we) — 
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whence 


Then inserting in (20) and (21) the estimates m. and s, of ws and a: 
calculated from the control series, and the first and second moments 
about the origin, m{, and mi, , of the treated series, we have 


(mis — ms) + + AL) 
PL) 


m, = Ms 


(221) 


(mig mism,) + (m2 (223i) 


2 
= & 


+ 


where 


L = (m. — m)/V — + 83 


If, then, we insert provisional estimates m, and s, into the righthand 
side of (22i), we get a closer approximation to the moment estimate of 
wu, , and using this, with the provisional s, , in (2211) we get a revised s, . 
The process may be repeated until stable values of m, and s, are reached. 
The proof of the convergence of this process is straightforward, and 
depends on an important property of the function A(), namely that 


0 < Xn) <1 


for all finite 7. This result is proved elsewhere (Sampford, 1952). 


6.4. The assumption of zero correlation in the ‘natural death’ model 


I have already suggested (§ 4.3) that when one of the two distribu- 
tions, F,(t,), results from deaths due to natural causes, and a control 
series is available from which an adequate estimate of the parameters 
of F,(t,) can be made, there is little point in attempting to estimate the 
marginal distribution F,(¢,), as this distribution could only operate in 
an idealised population with no natural death-rate, but that we should 
rather attempt to estimate a distribution G(t,) such that 


I(t, & gt) f(t). (23) 


The distribution G(‘) thus represents the additional death-rate 
superimposed on the natural death-rate by the action of the treatment. 
When /(é, , ¢,) is bivariate normal, it seems reasonable to try to approxi- 
mate to it by choosing G(é) as a normal distribution—in other words, 
to attempt to approximate by another bivariate normal distribution 
with the same values p, , ¢, , but having p = 0, regardless of the true 
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value of pin (4, , &). The justification for such a procedure can only 
be empirical; when the joint action is independent the results appear 
to be satisfactory. For example, from the data of Table V, 


m, = 10.07, & = 3.16; mi, = 9.06, mi, = 87.2340. 


of , have been calculated, by the method of the previous 
section, for a number of different values of p. The curves of the joint 
distribution functions, transformed to N.E.D.’s, are plotted in Figure 


1 for selected values of p; the curves corresponding to p = —0.75 and 
p = 0.75 are barely distinguishable, and those for larger values of 


pare not very different, suggesting that an adequate fit could be obtained 
by assuming p = 0, whatever its true value. This conclusion is supported 
by the results given in Table IV. The values uw, and o, naturally differ 
very considerably; however, the third and fourth moments of the joint 
distribution, us, and py, , show very little variation, and are practically 
constant over the range — 0.5 < p < 0.9. The quantities y, = us) /o7 
and y2 = (u4,/o3) — 3 are also tabulated. 


TABLE IV 
VALUES OF py and o,, AND CORRESPONDING VALUES OF psy , ayy . AND RELATED 
FUNCTIONS, ESTIMATED FROM DATA OF TABLE V, ASSUMING A RANGE OF VALUES 


OF p. 

p | Mas | x2 
—1.00 | 14.7 3.69 —9. 881 —0.845 89.4399 | 0.374 
—0.75 | 13.05 2.69 —8.258 | —0.706 | 92.8440 | 0.502 
-0.50 12.08 2.15 | —7.688 | —0.658 | 94.4717 | 0.562 
—~0.25 | 11.46 1.81 | —7.453 | —0.638 95.1872 | 0.590 

0 10.99 1.66 —7.362 | —0.630 | 95.5678 | 0.603 

0.25 | 10.62 1.55 | —7.329 | —0.627 95.6930 | 0.608 

0.50 | 10.24 1.52 | —7.328 | —0.627 | 95.7732 | 0.610 

0.75 | 9.84 1.58 | —7.300 | —0.624 | 95.9010! 0.616 

0.90 | 9.54 | 1.71 | —7.172 | —0.613 96.5374 | 0.639 

1.00} 9.23 | 1.96 | -6.425 | —0.550 | 100.2623! 0.779 


| 
| 


The corresponding quantities calculated from the sample, with ap- 
proximate standard errors, are 


= —7.5388 + 5.938 = —0.645 + 0.389 
my, = 88.5712 + 19.5751 y = 0.339 + 1.546 


No general formulae for variances of y, and y, are available, but they 
will presumably be of much the same order as those for a normal parent. 
For a large n these are, approximately, 
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Van) = 6/n; V(y2) = 24/n 


values which, taken in conjunction with the values of Table 1V, suggest 
that a very large sample would be required before a non-zero value of 
p would produce a joint distribution with skewness and kurtosis sig- 
nificantly different from those for the corresponding uncorrelated dis- 
tribution with the same mean and variance. 

If the joint action were dependent the bivariate distribution of 
t, and ¢, , even if normal, would probably not be such as to have F(t.) 
as marginal distribution. The ‘superimposed’ distribution G(¢) of (21) 
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Distribution functions of survival-time for the independent correlated joint action of two causes 
of death. 

In each case the two survival-times are assumed to have a bivariate normal distribution with 
correlation p: the mean and variance of one survival-time are kept constant throughout; those for the 
other are chosen to give a constant mean and variance, independent of p, to the resultant distribution 
function. (See §§ 6.2-6.4 and Table IV). 
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could be interpreted as the distribution of survival-times for all animals 
which die as a result of the treatment, or whose deaths from natural 
causes are accelerated by the treatment; the satisfactory fitting of the 
data by the assumption of a normal form for G(¢) seems less likely, how- 
ever, than for independent action. A plausible alternative form for 
G(t), when an assumption of normality proves inadequate, is desirable, 
but could only be based on the examination of a considerable number of 
samples sufficiently large to indicate departure from normality, and 
such data are not at present readily available. 

In gencral, then, it seems reasonable to suggest that, where natural 
mortality gives rise to survival-times normally distributed in some 
time-scale t, and a treatment is applied which might ideally be expected 
to give rise to a distribution of survival-times normal in the same time- 
scale, an attempt should be made to fit the resultant two-stimulus 
distribution by the form (23), where F,(t) is the distribution in the 
control series, and G(é) is some other normal distribution, and that this 
attempt may be made with some confidence if the joint action is in- 
dependent, and with at least a hope of success in the case of dependent 
action. It would, of course, be mathematically possible to assume the 
form (18), and to estimate p as well as ny, and o,. This procedure might 
lead to a rather better fit than the assumption of (23); the variances of 
the estimated parameters, however, would probably be large; if the 
joint action were dependent no very precise meaning could be attached 
to the estimated p; and even if the joint interaction were independent 
the estimates of p, u, and o, would probably be. of less practical use 
than the parameters of the ‘superimposed’ distribution G(¢), provided 
its assumed form gave an adequate fit. 


6.5. Estimation assuming zero correlation. 


We assume that the distribution of survival times in the treated 
population is 


Qn =1- H(t) 


where 


Q; = | du 


V 


and (1 — Q.) = F,(0) is the distribution function of survival-times in 
the untreated population. The data consist of two samples of survival- 
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times, one control series, from which estimates of (wu. , o2) can be made, 
and one treated series, from which we hope to estimate (u, , ¢,). The 
estimation proceeds in three stages: 


(i) Calculation of provisional estimates. 


For a few suitable values of ¢, Qy is observed directly from the 
treated sample, or read from a smooth curve drawn through selected 
sample points, and Q, either observed from the control series or ealeu- 
lated from the estimates of , and o. Q, is then calculated as the ratio 
of Q, and Q, , transformed to N.E.D’s, and plotted against ¢. A straight 
line is drawn by eye to fit the points; this line, 


=at+ byt 


is a provisional estimate of 
m =a, + Bit; 


from it the values a, and 6b, are determined, and provisional estimates 
of w, and o, are calculated as 


(m), = —a,/b, ; = 1/b, 
(ii) Estimation by moments. 


The iterative procedure of § 6.3, with p set equal to zero, is used 
to revise the provisional estimates, giving rise to value (7m,), ; (S:)> . 
(Alternatively stage (i) may be omitted, and values of m, and s, guessed, 
to provide a starting point for the iterative process; this is not recom- 
mended, however, as a bad guess may increase the time and labour 
required to carry through the iteration by far more than the amounts 
saved by omitting stage (i)). 

The estimates m, , s; , m, , 8 available at this stage are of varying 
efficiencies. The information on yz, and o, contained in the treated series 
has been ignored, so that m, and s, are not of maximum efficiency; 
however, all the information from the control sample has been used, 
and, provided the treated sample is not much larger than the control 
sample, the loss of efficiency will not be serious. More important, 
however, are the estimates m, and s, ; no general rule for their efficiency 
can be given, as a general form for the minimum variance cannot be 
evaluated. Clearly, where u, is very much less than p, , so that nearly 
all individuals in the treated series are killed by the treatment, the 
moment estimates m, , s, will have high efficiency, which will decrease 
as mw, increases towards p . Particular cases appear to suggest that 
when the two means are about equal, m, is moderately efficient, but 
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s, may have very low efficiency in the examp!ec of § 6.6 (in which m, 
is actually greater than m.) the variance of the moment estimate of 
o, is about ten times that of the maximum likelihood estimate. This 
surprisingly low efficiency appears to be due to the fact that o, appears 
in 19 (i) and (ii) only in the forms (07 + wi) and Vo; — 2po., + 03, 
both of which are rather insensitive to variation in ¢, when, as happens 
in this case, both w, and ¢, are of appreciably greater magnitude than o, . 

In view of this it seems safe to assert, as a general principle, that the 
moment estimates should not be accepted as final, but used as provisional 
estimates in the maximum likelihood method of the next section. For 
very large samples, however, and in rather more favourable cases than 
the example considered, the moment estimates may be sufficiently 
precise to warrant their being accepted as final, thus avoiding the 
necessity of performing the somewhat laborious maximum likelihood 
calculations (note, however, the qualifying remarks of § 6.10). In such 
a case, variances and covariances of the moment estimates can be 
calculated as follows. The values m. and s, will have been caleulated 
as the mean and standard deviation of the control series; they will thus 
be uncorrelated, and have error variances estimated by 


Vom.) = s3/n; V(s,) = s3,/2n 


where n is the size of the control sample. The remaining variances and 
covariances may be estimated approximately (for large samples) by 
deriving, from the two moment formulae used in the iterative procedure, 
estimates of increments of m, and s, in terms of increments in the four 
observables m, , s, , mj, and mj, . and calculating expected values of 
squares and products of increments. The resultant formulae are: 


A’? V(m,) = BiV(m',) + 2B,Be cov (mi, , mz) 
+ + BEV(m.) + (s) (24i) 

s,4° cov (m, , = B,C, + (B,C. + BC,) eov (mi, , 
+ BC LV + Om.) + BYC,V(s,) (2411) 

si V(s,) = CTV (mis) + cov (mis, mis) 
+ + CEVOn,) + (24iii) 


A cov (m,, m.) = B,V(m.) A cov (m, , = 
(2-4iv) 


cov (s; ms) 


| 
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where 
A = P(1 — siy) 
B, = (1 -— C, = —(m, — s¢) 
B, = 3¢ C,= 3 
= v) C, = (Ps; — — Q(m, — 
Ps; — 
B, = (P — Q¢s: Cy = 
La 
V (si + 82) 
(i.e. the value calculated in the last cycle of the iteration procedure) 
P,Q = P(D), QL) 
(s; + 82) V (si + 82) 
+ 
and Y= 


The variances and covariance of mj, and m3, , which appear in these 
formulae, are given by 


V(mis) = (uss — 


1 
cov , m5) = (uss — 


n 
Vim.) 
where n’ is the size of the treated sample. In these formulae u{, and 
uz, may be replaced by m{, and m3, , while us, and ui, may either be 
estimated directly by the third and fourth sample moments about the 
origin, or calculated from formulae (19) as functions of m, , s; , M2, 
and s, . 

In view of the remarks above about the possible inefficiency of the 
moment method, it may be suggested that it is a waste of time, and 
that a more profitable procedure would be to use the initial approxima- 
tions, a, and b, of (i) above, directly in the maximum likelihood process, 


- 
wit 
a 
si 
si 
| 
jane 0 
si 
st 
fi 
t 
i 
€ 
Pie 
; 
Wi 
- 
| 
| 
ft 
ie 
i 


RESPONSE-TIMEE DISTRIBUTIONS 337 


No general ruling can be given on this point: each case must be con- 
sidered on its own merits, taking into account previous experience of 
similar analyses, and the relative magnitudes of the provisional estimates 
of the parameters. However, the moment procedure is unlikely to 
seriously worsen the approximations (even in the particularly un- 
favourable example here given the estimates a, and b, calculated from 
the moment estimates are rather closer to the final values than the 
initial approximations), and it may well improve them sufficiently to 
reduce the number of cycles required in the maximum likelihood process 
—a saving which, even for a reduction of only one cycle, will more than 
outweigh the time spent on the moment estimation. 


(iii) The maximum likelihood procedure. 


For simplicity of the resultant formulae this process is developed 


in terms of a = — p/o and B = Iya. We have 
= 
h(t) = +- 8.7 
where 


and we shall denote v(n,), by , A, Where = 1,2. Then 
log L => log (8;7,Q2 + 


where L is the likelihood of the treated sample only, and summation is 
over survival-times of all animals in the treated sample. Then 


Ada, 1Q2 + 


log L — BimtZ,Qo — B2tZ, 
Op, 


da, + 8 7 


Alog _ { 
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(This Z should not be confused with the Z, above.) 
As in § 5.2 the equations 


log L 


cannot be solved explicitly, so we have recourse to an iterative process 
involving four equations similar to equations (8). Generally 


log L _ (2 log log 1 
0009 00 dg + Li 

and the second term may be replaced by its expected value, zero. The 
expected value of the first term cannot be evaluated, but it may be 
calculated from provisional estimates, and used to improve the ap- 


proximation in the usual way. The maximum likelihood equations 
then become 


>, W? + 5b, WX + ba, D WY + 5b. WZ 
ba, WX + 8b, DX? + Ga, + 8b, XZ 
da, WY + 6b, D XY + 8a, Y? + th YZ =-DY 
da, WZ + 6b, XZ + ba, + bb. DZ -> 


where W, X, Y, and Z are calculated from first approximations a, , b, , 
a, , b. , and da, , 5b, , da, , 5b, are the corrections to be applied to these 
approximations. These equations are based only on the treated sample, 
and are appropriate if no control series is available. In general, how- 
ever, a control series will have been run, and this adds to the log likeli- 
hood a term 


x 


log L, = n log — (ar + BoA)? 


where x is the size of the control scries, and }~, denotes summation 
over survival-times of all members of this series. Then 


2 


log L, log L, ® log L, 


giving, for the maximum likelihood equations based on both samples 
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sa, + oh, DWN 4+ fa, 
+ 6, DwWZ=-DW | 
+ 6h, 
da, >, WY + 6b, DENY 4+ fn + 
+ dh, > ¥Z} = —(na, + 


(25) 


da, WZ + 6b, XZ + YZ! 


+ + + > = -(a, Pt+b, 


2 
2 


The calculations are most conveniently performed as shown in 
Table VIII. ¢is tabulated, (together with n/ , the number of deaths in 
the treated series at time ¢, if mj is not always 1), and y, , 2 . estimates 
of m , mn , enleulated from the provisional estimates. v, and vy, are read 
from tables, together with \, . A. , which are most conveniently tabulated 
direct as B,A,; Finally (b,», + is tabulated, and Y. and 
Z and their sums and sums of squares and products calculated. 

The equations are most conveniently solved by inverting the matrix 
of coefficients; if the adjustments 6a, ete. are large the process may be 
repeated until stable values are attained. The inverse of the coefficient 
matrix for the fihal cycle provides estimates of the variances and covari- 
ances of the estimated parameters. 


6.6. Example. 


The data of Table V are survival-times, in units of 2 months, of 
two groups of mice; one control series and one series subjected to con- 
tinuous whole-body gamma-radiation at an intensity of 2.2 standard 
units of radiation daily. The data are abstracted from a paper by 
Lorentz ef al. (1917) and grouped for the purposes of the example. On 
this scale, treatment actually began at time 7 = —1.5. 
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TABLE V 


SURVIVAL-TIMES (1 UNIT = 2 MONTHS) OF MICE IN ONE CONTROL SERIES AND 
ONE SERIES TREATED WITH GAMMA-RADIATION (2.2 UNITS DAILY). 


Survival time | Numbers surviving for time 7 in 


Control series | Treated series 


| 


= 


The cumulative mortality curves, transformed to N.E.D.’s, are plotted 
against time in Figure 2; from the apparent linearity of the ‘control’ 
line, the similarity of the ‘treated’ curve to the curves of Figure 1, and 
consideration of other series of data from the same source, it appears 
that survival-time in the control series is normally distributed, and 
that an assumption of a normal form in untransformed time for the 
‘superimposed’ distribution should give an adequate fit. 
From the control series 


n=58 584 = 6450 


giving provisional estimates 
m, = 10.07 8 = 3.134 


(As the estimate s, is to be improved by a maximum likelihood 
procedure, it has been calculated as — this, the maximum 
likelihood estimator from the control series only, seems likely to approxi- 
mate more closely to the final value than the unbiased estimate 
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1). The latter is, of course, preferable if only the first 
tivo stuges of estimation are to be used, so that s, so calculated is the 
fina! estimate of 

From the treated series 


n = 47 Woe 9.06 
= 42 87.2340 
Lf 4100) mi, = 876.127660 


= 41178 mi, = 9089.78723404. 


= 427220 
TREATED 
CONTROL 
x 
x 
x 
NED 
of 
% 
Kill 
Of 
+ 
“15+ 
+ 
x 
Time Ciont = 2 moth) 
FIG. 2. 


Mortality among mice in a control series and a series subjected to gumma-radiation. 
Data from Lorenz et ol., 1947. Points represent observed mortality; the curves are calculated 
from the parameters estimated in the example of § 6.6. 
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Using the estimates m, and s, to calculate y. , and reading yy, from a 
smooth curve drawn by eye to fit the points of Figure 2, the provisional 
estimation is completed as in Table VI 


TABLE VI 
CALCULATION OF PROVISIONAL E“TIMATES FROM DATA OF TABLE V 
t Y2 UH Qu 
7.5 —0.82 0.7939 —0.8 0.788 0.993 —2.46 
8.5 —0.50 0.6915 —0.4 0.655 0.947 —1.62 
9.5 —0.18 0.5714 0.1 0.460 0.805 —0.86 
10.5 0.14 0.4443 0.6 0.273 0.614 —0.29 
11.5 0.46 0.3228 i 0.115 0.356 0.37 
12.5 0.7 0.2177 =| 1.9 0.029 0.133 
| | 


The resultant values of y, give a fairly good straight line when plotted 
against time, and a line drawn by eye gave the provisional estimates 


a, = —7.63; b, = 0.70, 
and hence 
m, = 10.90 s, = 143 


These estimates of m, and s, are now improved by the moment 
method. Assuming p = 0, and inserting into formulae (22) the values 
given above for m, , s. , m{, and m3, , the iteration equations become 


Vs; + 9.8220 Z(L) — 1.01 


m, = 10.07 + P(L) 
1.01m, — 13.8822 
s; = 9.8220 + — 
where 
_— 10.07 — m, 


V(si + 9.8220) 
The details of the iteration are shown in Table VII, and give rise to 
estimates 
m, = 10.98 s, = 1.672 
which, with 


Mm, = 10.07 Sg 3.134 
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TABLE VII 


my V si +9 . 8220 L | Z 
10.90 2.0419 | 3.445 —0.24 0.3876 0.4052 
10.87 —0.23 0.4090 
2.7230 | 3.542 —0.23 0.3885 0.4090 
10.97 0.4013 
2.8381 | 3.558 | —0.25 0.3867 0.4013 
10.98 : —0.26 0.3974 
2.7953 | 3.552 0.3857 0.3974 
10.98 | —0.26 0.3974 
2.7953 | 
| 
| 
i 


are used to give provisional estimates for the maximum likelihood 
method: 


a, = — 6.57, by = 0.598; a = —3.21, be = 0.319. 


Had it been intended to stop the analysis at this point, s, would 
have been chosen, as suggested above, as the unbiased estimate 3.162, 
which would have produced slight differences in m, and s, . However, 
for the sake of the example, it is convenient to calculate error variances 
and covariances for the present estimates by the methods of the previous 
section. We have 


V(m,) = 0.1693; cov (Mz , $2) = 0; V(s.) = 0.084672 


L = —0.26 »(L) = 0.6401 ML) = 0.5761 
v(—L) = 0.9704 \(—L) = 0.6894 
whence formulae (24) give 
V(m,) = 0.3749; cov (m, , &) = —0.03304; V(s,) = 3.577779 
cov (m, , m2) = —0.0969; cov (m, , s:) = —O.044118. 
cov , M2) = 0.34719; cov , = —0.255432. 
Finally, the provisional estimates given above may be improved by 


the maximum likelihood method; details of the analysis are shown in 
Table VIII. The equations for the correction terms are: 
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9.7176 6a,4+- 103.2105 6b,+ 10.5562 6a,4+ 124.9590 6b, = --0.231 
103.2105 6a, +1116.8413 66,4 102.9067 64,4 1253.5386 6b, = —1.596 
10.5562 6a,+ 102.9067 6b,+ 91.2391 6a,4+ 858.3505 6b. = —0.120 
124.9590 da,+ 1253.5386 6b,4+858.3505 6a,+9553.6588 6h, = 1.262 


The inverse of the coeflicient matrix is 
6.01963698 —0.54967381 —0.09233504 0.00168368 | 


—0.54967381 0.05125006 0.00915917 —0.00035787 


(V] = (26) 
—0.09283504  0.00915917 —0.07275659 —0.00652910 


0.00168368 —0.00035787 —0.00652910 0.00071621. 


and 
| r—0.231] [—0.50430) 
5b, —1.596 0.04453 
ba, ~0.120} | 0.00622 
| 5b. J |—1.262} | 0.00006] 
whence 
a, = —7.07 +245 a, = —3.20 +0.27 
b, = 0.643 40.226 0.319 + 0.027 


The corrections to a, and b, amount only to some 20% of their 
standard errors, while a, and 6, are practically unchanged. Clearly 
another cycle is unnecessary. The test of goodness of fit is shown in 
Table IX. The x’ test should strictly be applied simultancously to the 
control and treated series, the two x” values being added and tested with 
degrees of freedom 6 fewer than the total number of cells in both series 
(for the deviations in each series sum to zero, so that. 1 degree of freedom 
must be deducted for cach series, and four parameters have been esti- 
mated from the sample). The treated sample alone has provided all the 
information on two parameters, and a part of the information on two 
more; the x’ calculated from this series alone may therefore be supposed 
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TABLE IX 
TEST OF GOODNESS OF FIT ON DATA OF TABLE V. 


t “ y 1 Q2 On Exp. no. Obs. no. x? 
killed killed 
3.5 | —4.82] -—2.08 | 1.0000 | 0.9812 | 0.9812 | 
4.5 | —4.18 | -—1.76 | 1.0000 | 0.9608 | 0.9608 | 0.961. . 1 
5.5 | -3.53! -1.45 | 0.9998 | 0.9265 | 0.9203 | 1.67/52 2/° 
6.5 | -2.39] -1.13 | 0.9981 | 0.s708 | 0.8601 | 2.69) 2) 
7.5 | —2.25 | -—0.81 | 0.9878 | 0.7910 | 0.7813 | 4.13 4 0.004 
8.5 | -1.60] -0.49 | 0.9452 | 0.6879 | 0.6502 | 6.17 7 0.112 
9.5 | —0.96 | -—0.17 | 0.8315 | 0.5675 | 0.4719 | 8.38 8 0.017 
10.5 | —0.32 0.15 | 0.6255 | 0.4404 | 0.2755 | 9.23 10 0.064 
11.5 | 0.32 0.47 | 0.3745 | 0.3192 | 0.1195 | 7.33 4 1.513 
12.5 | 0.97 0.79 | 0.1660 | 0.2148 | 0.0356 $041 5 ot ts 1.018 
! 
| | Total 2.93 


to have degrees of freedom greater than (i: — 5) and less than (/ — 3), 
where k is the number of cells in the treated series. Thus if this x’ is 
not significantly large when tested with (k — 5) degrees of freedom, the 
normal form assumed for G(/) may be safely accepted; if it is significantly 
large when tested with (k — 3) degrees of freedom, the assumed (Gt) 
must be rejected (always provided the assumption of normality for 
F(t) is justified). If the observed x’ lies between the significance levels 
of (k — 5) and (k — 3) degrees of freedom, the complete test must be 
made. In the present example, the treated series alone gives a x” (based 
on 7 cells) of 2.73; tested against 2(= 7 — 5) degrees of freedom, this is 
not significant; apparently the assumption of normality for G(t) gives 
an adequate fit. For the complete test it is necessary to combine this 
x’ with that for the control series; this value, for which details of calcula- 
tion are not shown, is 4.89, based on 9 cells, giving a total x’ of 7.62, 
with 10(= 7 + 9 — 6) degrees of freedom. 

The cumulative curves for the fitted distributions, transformed to 
N.E.D.’s, are shown in Figure 2. 

Approximate estimates of the variances of the maximum likelihood 
estimates of the m’s and s’s can be calculated from those of the a’s and 
b’s for comparison with those of the moment estimates. For 


= 1/b; 


and 
m, = —a,/b, 


so that, approximately, 
V(s,) = V(b,)/b: 


a 
t= 
; 
; 
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and 
Vom) = {V(a) + 2m, cov (a; , b,) + V(b 
From these formulae 
(m,) = 0.3098; Vis,) = 0.299812; V(m,) = 0.1359; V(s,) = 0.069163 


Comparison of these values with those given above for the moment 
estimates shows that while the moment estimates for m, , m2 , and s, 
have quite high efficiencies, that for s, is very inefficient. This point 
has already been discussed in § 6.5. 


6.7. A possible simplification of the maximum likelihood procedure. 


If the control sample is large, so that the moment estimates m, and 
s, , and hence a, and b, , are fairly reliable, some of the labour of the 
maximum likelihood procedure described above may be eliminated by 
accepting as final the values a, and b, claculated from the control series, 
and maximising the likelihood of the treated sample for a, and £, only, 
giving equations. 


sa, +54, DWX=-DW 
da, >, WX 4+- 5b, = 


where W and X are the functions defined in the previous section. This 
has the disadvantage that the variances of a, and 6, given by the inverse 
matrix of coefficients are biased, being smaller than the true maximum 
likelihood variance estimates, and that no estimates are available of 
covariances between (a, , b,) on the one hand and (a, , b2) on the other. 
Applied to the example of the previous section, this method gives 
estimates 


II 


—7.03 + 2.36 (ly 
b, = 0.640 + 0.220 bs 


—3.21 + 0.33 
0.319 + 0.030. 


Comparing these with the results of the previous section, the bias in 
the standard error of a, and b, is seen to be negligible. The standard 
errors of a, and b, (calculated by asymptotic formulae from V(m,) and 
V(s.)) are larger than those of the previous section; this is to be ex- 
pected, as the information on these parameters contained in the treated 
series has been neglected. 


a; 


Il 


6.8. Estimation from grouped data. 


The problems oceurring in the estimation of parameters from 
grouped response-time data have been discussed in some detail in the 
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first paper of this series. In the present case there are two alternative 
approaches. If the grouping interval is small (say, less than the smaller 
of o, and o,) all deaths in any interval may be considered as having 
occurred at a class-mark in the interval, and the method for ungrouped 
data used (as has, in fact, been done in the example of § 6.6). If the 
grouping interval is larger a special maximum likelihood procedure 
taking account of grouping may be developed, similar to that given in 
the previous paper for the normal distribution (Sampford, 1952; § 3.7). 
In practice the method of § 6.5, applied to class-marks, may conveniently 
be used in all cases up to the estimation of », and o, by moments; at 
this stage provisional estimates of both o, and o, are available, and the 
final analysis may be made by the ‘ungrouped’ or the ‘grouped’ maximum 
likelihood method according as these estimates are both greater than 
the grouping interval or not. 


6.9. Estimation of functions of the parameters. 


If a, , b, , a, , b, are estimates of the a’s and £’s, the estimated pro- 
portion killed at time 7, when both causes of death operate, is given by 


pr (27) 


Exact expressions for the variance and fiducial limits of pry are not 
available but its variance may be estimated, approximately, as 


V(pr) = (1 — pr)*J(T) (28) 
where 
J(T) = V(a,) + 27 cov (a, , + T°V(b,)} 
+ 2v,v.{cov (a, , a2) + T cov (a, , 
+ T cov (az , b;) + T’ cov (b, , be)} 
+ 3{V(az) + 27 cov (az , + T°V(b.)} (29) 
where 
= v(a; + b,T). 
The time fp. , required to kill a proportion P’ is given as the solution 
of 
Q(a, + + betp-) = 1 — P’ = Q’. (30) 
This may most conveniently be solved by obtaining a provisional 


estimate from the cumulative sample N.E.D. curve, and calculating 
successive increments by the formula 


bi, + (31) 
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where the Q, , v, are calculated from a provisional estimate. ‘The variance 
of dé; , may be estimated, approximately, as 


where J(t) is the function of (29), and»; = v(a; + b,tp-). Consider, 
for example, the data of § 6.6. The estimated parameters are 

a, = —7.07, b, = 0.648, a, = —3.20, b, = 0.319, 
and their variances and covariances are given as elements of the matrix 


{[V] (formula (26)). The proportion killed at 18 months is given by 
(25) with 7’ = 7.5. 


a, + b,7 = —2.25 Q, = 0.988  », = 0.032 
a2 +bT = —081 Q,=0.791 » = 0.363. 
Then 
Pr =1—- QQ, 
= 0.218 
J(7.5) = (0.001024 X 0.657346) + (0.023232 x —0.031144) 
+ (0.131769 X 0.015087) 
= 0.00193758 
S.E.(p;) = 0.782-VJ(7.5) 
= (0.034. 


Again, the time required to reach a 50% kill is given by Q,Q. = 0.50. 
Taking a provisional é.; at 9.0, successive approximations are 9.4, 
9.36, and, finally, 

tp, = 9.35 (i.e. 21.7 months after the beginning of treatment) 
Then 


a, + btp- = —1.06; vy, = 0.266; 
a, + betp, = —0.22; v, = 0.663. 
(0.070756 X 0.221145) + (0.352716 X —0.022240) 

+ (0.439569 X 0.013255) 


J (lp) 


II 


0.01362942. 


x 0.1167 
= 3gt7 


).31. 


i 
ote 
q 
j 
tied 
i 
f 
i 
q 


350 BIOMETRICS, DECEMBER 1952 


These standard errors are useful as giving a rough idea of the pre- 
cision of the estimates py and fp. ; that for tp- could be used to attach 
limits to the estimate, although, in the absence of any knowledge of 
the form of the distribution of /». , such limits can only be of doubtful 
validity. V(pr) should never be used to calculate limits, as p, will 
probably be far from normally distributed in nearly all cases. 


6.10. Concluding remarks on estimation with control mortality. 


The analysis of data of this type has been developed at some length 
in the preceding sections, and the experimenter who wishes or is obliged 
to perform his own statistical analysis may well feel that the amount of 
computing required is excessive. However, this approach is only 
necessary where there is a considerable mortality among the controls, 
and this is likely to occur only in comparatively long-term experiments. 
In consequence the analysis, though admittedly tedious compared 
with, say, that of a routine assay, will still account for only a small 
fraction of the time and money spent on the total investigation. For 
example, the data of Table V consist of observations made at regular 
intervals over a period of nearly three years; compared with this the 
time required to carry through the calculations of § 6.6 is negligible. 


6.11. Two applied treatments. 


If the two causes of death are two applied treatments, rather than 
one treatment and some unavoidabie cause resulting in what I have 
called a ‘natural’ death-rate, the situation is greatly simplified in sc 
far as it is possible to estimate the parameters of the two distributions 
of survival-times by applying the treatments separately. They are 
thus likely to be applied together only in an attempt to investigate the 
nature of their joint action. The possible types of joint action have 
already been discussed in §§ 4.1-4.3. Similar action is only likely to 
occur when the two treatments are doses of poisons having identical 
or similar active constituents, or at least having a common site of 
action. The consequent distribution of survival-times is that due to 
the effective combined dose, and so might be expected to be normal in 
terms of that time-metameter which normalizes the two separate dis- 
tributions. Dissimilar action, on the other hand, would produce a 
non-normal joint distribution; thus a test of normality should provide 
a test for similar as opposed to dissimilar action. 

On the hypothesis of independent action with zero correlation, the 
mean of the joint distribution is given (§ 6.2) by 


wis = mP(A) + — Voi + 03 Z(A) 
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with 


For correlation p > 0, the mean of the joint distribution is greater than 
this quantity, while for p < 0, the mean is less. This suggests that a 
suitable test for zero correlation might be based on the quantity 


= mi, — mP(L) — mQ(L) + Vs, + 8; ZL) 


where m, , , 8 , 8 are estimates of , , , calculated from 
series in which the treatments are applied individually, and L is obtained 
by substituting these estimates in A. On the hypothesis p = 0, é has 
for large samples a zero expected value and an asymptotic variance 
estimated by 


2 2.2 2 2 4 4 
= 4 4 (2 +2) 


No + \n, ne 


where sj; is the variance of the joint distribution, P, Q, and Z are func- 
tions of L, and n, , n. , n, are the sizes of the three tested samples. 
Further ¢ may probably be assumed normal without introducing too 
serious an error. A significantly large positive would then indicate 
a positive correlation, and a negative ~ a negative correlation. If an 
actual value of p were required, it could then be estimated by using 
the form for yu{,; given in § 6.2, inserting the estimates mj, , m, , 8, , m2 
and s, , and solving the resultant equation in p by trial and error or by 
an iterative process. This procedure, of course, ignores the possibility 
of dependent action; once this possibility is admitted the situation 
becomes extremely complicated, and there does not seem to be any 
method, depending only on survival-times, that will distinguish all the 
possible forms of dependent action. When the survival-times are 
uncorrelated, but each treatment accelerates the effect of the other, a 
negative ~ should result; if each treatment retards the action of the 
other, should be positive. The further complications introduced when 
the effect of one treatment accelerates, but is itself retarded by, the 
action of the other, or when the joint action is dependent with corre- 
lated survival-times, are beyond the scope of this test. However, it 
is just conceivable that for certain pairs of treatments, prior knowledge 
might enable the experimenter to say definitely that if the response- 
time were correlated, p was positive, and that the only plausible type of 
interaction was mutual acceleration. In this case a positive ¢ would 
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imply a positive correlation of survival-times sufficient to outweigh any 
mutual acceleration, while a negative — would indicate a mutual ac- 
celeration sufficient to overcome the effect of any correlation of survival- 
times. A value of £ not significantly different from zero might imply 
that the joint action was independent and uncorrelated, but might 
also mean that the increase in mj, produced by a positive correlation. 
and the decrease produced by mutual acceleration of the treatment 
effects, were of the same order of magnitude: in short, the test based on 
t, though incapable of discriminating between all acceptable types of 
joint action, might in this case serve to detect one important feature, 
without necessarily ruling out the other. It is clear that this approach 
is far from perfect, but little improvement seems possible, in view of 
the unpleasantly complicated nature of the distributions involved, and 
the difficulties of producing a practicable mathematical model for 
dependent action and of distinguishing, by the form of the response- 
time distribution alone, between uncorrelated independent. action, 
correlated independent action, and dependent action. 


SUMMARY 


The types of response time distribution occurring when two or more 
stimuli act on a sample of individuals are discussed. Two situations 
are discussed in considerable detail; the ‘accidental death’ model, in 
which each response can be related to its appropriate stimulus, and the 
‘natural death’ model, in which the exact cause of any death cannot be 
determined, but the distribution of ‘potential survival-times’ to the 
two stimuli can be assumed bivariate normal. Maximum likelihood 
methods of estimation are developed for these situations, and tables 
are given to simplify the calculations. 

I am indebted to Dr. D. J. Finney, Dr. H. O. Hartley, and Dr. I. 
Sutherland for helpful criticism, to Dr. R. M. Mole, who first drew my 
attention to the problems involved in the estimation of multi-stimulus 
response-time distributions, to Dr. D. B. Carlisle and Dr. P. F. R. 
Dohrn for permission to use the data of Table I, and to Miss Mary 
Callow for her assistance in preparing the Appendix Tables. 
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APPENDIX 


Tables of the functions », A, and ¢. 


| | 
n v | 
--5.0 0.00000 0.00001 0.00004 
—4.9 0.00000 0.00001 0.00006 
-4.8 0.00000 0.00002 0.00010 
-4.7 0.00001 0.00003 0.00015 
—4.6 0.00001 0.00005 0.00022 
—4.5 0.00002 0.00007 0.00034 
—4.4 0.00002 0.00011 0.00051 
—4.3 0.00004 0.00017 0.00075 
—4.2 0.00006 0.00025 0.00110 
-4.1 9.00009 0.00037 0.00159 
—4.0 0.00013 0.00054 0.00227 
-3.9 0.00020 0.00077 0.00322 
—3.8 0.00029 0.00111 0.00451 
0.00042 0.00157 0.00624 
—3.6 0.00061 0.00220 0.00855 
—3.5 0.00087 0.00306 ! 0.01157 
—3.4 0.00123 0.00419 0.01549 
—3.3 0.00172 0.00569 0.02050 
—3.2 0.00239 0.00764 0.02683 
—3.1 0.00327 0.01015 0.03473 
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v 
—3.00 0.00444 0.01333 0.04444 
—2.99 0.00457 0.01369 0.04552 
—2.98 0.00471 0.01406 0.04662 
—2.97 0.00485 0.01444 0.04774 
—2.96 0.00500 0.01483 0.04889 
—2.95 0.00515 0.01522 0.05005 
—2.94 0.00531 0.01562 0.05124 
—2.93 0.00546 0.01604 0.05245 
—2.92 0.00563 0.01646 0.05369 
—2.91 0.00579 0.01689 0.05494 
—2.90 0.00596 0.01733 0.05622 
—2.89 0.00614 0.01778 0.05752 
—2.88 0.00632 0.01824 0.05885 
—2.87 0.00650 0.01871 0.06020 
—2.86 0.00669 0.01919 0.06157 
—2.85 0.00689 0.01968 0.06297 
—2.84 0.00709 0.02018 0.06439 
—2.83 0.00729 0.02069 0.06584 
—2.82 0.00750 0.02121 0.06731 
—2.81 0.00772 0.02174 0.06881 
—2.80 0.00794 0.02228 0.07033 
—2.79 0.00816 0.02284 0.07188 
—2.78 0.00839 0.02340 0.07345 
—2.77 0.00863 0.02398 0.07505 
—2.76 0.00887 0.02457 0.07667 
—2.75 0.00912 0.02517 —0.07832 
—2.74 0.00938 0.02578 —0.08000 
—2.73 0.00964 0.02640 —0.08171 
—2.72 0.00990 0.02704 —0.08344 
—2.71 0.01018 0.02768 —0.08520 
—2.70 0.01046 0.02834 —0.08699 
—2.69 0.01074 0.02902 —0.08880 
—2.68 0.01104 0.02970 —0.09064 
—2.67 0.01134 0.03040 —0.0925! 
—2.66 0.01165 0.03111 —0.09441 
—2.65 0.01196 0.03184 —0.09633 
—2.64 0.01228 0.03258 —0.09828 
—2.63 0.01261 0.03333 —0.10027 
—2.62 0.01295 0.03409 —0. 10228 

0.01329 0.03487 —0. 10432 
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n v 
—2.60 0.01365 | 0.03567 —0. 10638 
—2.59 | 0.01401 0.03648 —0.10848 
—2.58 0.01438 0.03730 —0.11060 
—2.57 0.01475 0.03813 —0.11276 
—2.56 0.01514 0.03898 —0.11494 
—2.55 0.01553 0.03985 —0.11715 
—2.54 0.01594 0.04073 —0.11939 
—2.53 0.01635 0.04163 —0.12166 
—2.52 0.01677 0.04254  —0. 12396 
—2.51 0.01720 0.04346 —0. 12629 
—2.50 0.01764 0.04441 —0. 12865 
—2.49 0.01809 0.04536 —0.13104 
—2.48 0.01855 0.04634 | —0.13346 
—2.47 0.01901 0.04732 —0.13591 
—2.46 0.01949 0.04833 —0.13838 
—2.45 0.01998 0.04935 —0.14089 
—2.44 0.02048 0.05039 —0.14342 
—2.43 0.02099 0.05144 —0.14599 
—2.42 0.02151 0.05251 —0.14858 
—2.41 0.02204 0.05360 —0.15121 
—2.40 0.02258 0.05470 —0.15386 
—2.39 0.02313 0.05582 —0.15654 
—2.38 0.02370 0.05696 —0.15926 
—2.37 0.02427 0.05811 —0.16200 
—2.36 0.02486 0.05928 —0.16477 
—2.35 0.02546 0.06047 —0. 16757 
—2.34 0.02607 0.06168 —0.17039 
—2.33 0.02669 0.06290 —0.17325 
—2.32 0.02733 0.06414 —0.17614 
—2.31 0.02797 0.06540 —0.17905 
—2.30 0.02863 0.06668 —0.18199 
—2.29 0.02931 0.06797 —0.18497 
—2.28 0.02999 0.06929 —0.18796 
—2.27 0.03069 0.07062 —0.19099 
—2.26 0.03141 0.07196 —0.19404 
—2.25 0.03213 0.07333 —0.19713 
—2.24 0.03287 0.07472 —0.20023 
—2.23 0.03363 0.07612 | 0.20337 
—2.22 0.03440 0.07754 20653 
—2.21 0.03518 0.07898 —0.20972 
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| 
n v | 

—2.20 0.03597 0.08044 | 0.21294 
—2.19 0.03679 0.08192 0.21618 
—2.18 0.03761 0.08341 —0.21945 
—2.17 0.03845 0.08493 —0. 22274 
—2.16 0.03931 0.08646 —0. 22606 
—2.15 0.04018 0.08801 —0.22941 
—2.14 0.04107 0.08958 —0.23278 
—2.13 0.04198 0.09117 —0.23617 
—2.12 0.04290 0.09278 —0. 23958 
—2.11 0.04383 0.09441 —0.24303 
—2.10 0.04478 0.09605 —0.24649 
—2.09 0.04575 0.09772 —0.24998 
—2.08 0.04674 0.09940 —0.25349 
—2.07 0.04774 0.10110 —0.25702 
—2.06 0.04876 0. 10282 —0.26058 
—2.05 0.04980 0. 10456 —0.26415 
—2.04 0.05085 0. 10632 —0.26775 
—2.03 0.05192 0.10810 —0.27137 
—2.02 0.05301 0.10990 —0.27501 
—2.01 0.05412 0.11171 —0.27867 
—2.00 0.05525 0.11355 —0. 28234 
—1.99 0.05639 0.11540 —0.28604 
—1.98 0.05756 0.11727 —0.28976 
—1.97 0.05874 0.11916 —0.29349 
—1.96 0.05994 0.12107 —0.29724 
—1.95 0.06116 0.12300 —0.30101 
—1.94 0.06240 0.12495 —0.30480 
—1.93 0.06336 0.12691 —0.30860 
—1.92 0.06494 0. 12890 —0.31242 
—1.91 0.06624 0.13090 —0.31625 
—1.90 0.06756 0.13292 —0.32010 
—1.89 0.06890 0.13496 —0.32397 
—1.88 0.07026 0.13702 —0.32784 
—1.87 0.07164 0.13909 —0.33173 
—1.86 0.07304 0.14118 —0.33564 
—1.85 0.07446 0.14329 —0.33955 
—1.84 0.07590 0.14542 —0.34348 
—1.83 0.07737 0.14757 —0.34742 
—1.82 0.07885 0.14973 —0.35137 

0.08036 0.15191 —0.35533 
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n v 
—1.80 0.08189 0.15411 —0.35930 
—1.79 0.08344 0. 15633 —0.36327 
-1.78 0.08502 0. 15856 —0.36726 
-1.77 0.08662 0. 16081 —0.37126 
—1.76 0.08824 0. 16308 —0.37526 
—1.75 0.08988 0.16536 —0.37926 
—1.74 0.09154 0.16766 —0.38328 
—1.73 0.09323 0.16998 —0.38730 
—1.72 0.09494 0.17232 —0.39132 
-1.71 0.09668 0.17466 —0.39535 
-1.70 0.09844 0.17703 —0.39939 
~1.69 0. 10022 0.17941 —0.40342 
—1.68 0.10202 0.18181 —0. 40746 
—1.67 0.10385 0.18422 —0.41151 
—1.66 0.10571 0.18665 —0. 41555 
—1.65 0.10759 0.18909 —0.41959 
—1.64 0. 10949 0.19155 —0. 42364 
—1.63 0.11142 0.19403 —0.42768 
—1.62 0.11337 0.19651 —0.43172 
—1.61 0.11535 0.19902 —0.43577 
—1.60 0.11735 0.20153 —0.43981 
—1.59 0.11938 0.20407 —0. 44384 
—1.58 0.12143 0.20661 —0. 44788 
—1.57 0.12351 0.20917 —0.45191 
—1.56 0. 12562 0.21174 —0. 45593 
—1.55 0.12775 0.21433 —0. 45995 
—1.54 0.12990 0.21693 —0.46397 
—1.53 0.13209 0.21954 —0.46798 
—1.52 0.13429 0.22216 —0.47198 
-1.51 0. 13653 0.22480 —0.47597 
—1.50 0.13879 0.22745 —0. 47996 
—1.49 0.14108 0.23011 —0. 48394 
—1.48 0.14339 0.23278 —0.48791 
—1.47 0. 14573 0.23547 —0. 49187 
—1.46 0.14810 0.23816 —0. 49582 
—1.45 0. 15050 0.24087 —0.49976 
—1.44 0. 15292 0.24359 —0.50368 
—1.43 0. 15537 0.24632 —0.50760 
—1.42 0.15785 0.24906 —0.51150 
—1.41 0. 16035 0.25180 —0.51539 
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” v | 
—1.40 0.16288 0.25456 —0.51927 
—1.39 0.16544 0.25733 —0.52313 
—1.38 0. 16803 0.26011 —0.52698 
—1.37 0.17064 0.26290 —0. 53082 
—1.36 0.17329 0.26570 —0.53463 
—1.35 0.17596 0.26850 —0.53844 
—1.34 0.17866 0.27132 —0.54222 
—1.33 0.18138 0.27414 —0.54599 
—1.32 0.18414 0.27697 —0.54974 
—1.31 0. 18692 0.27981 —0.55347 
—1.30 0.18974 0.28265 —0.55719 
—1.29 0.19258 0.28551 —0.56088 
—1.28 0.19545 0.28837 —0.56456 
—1.27 0.19834 0.29124 —0.56821 
—1.26 0.20127 0.29411 —0.57185 
—1.25 0.20423 0.29699 —0.57546 
—1.24 0.20721 0.29988 —0.57906 
—1.23 0.21022 0.30277 —0.58263 
—1.22 0.21327 0.30567 —0.58618 
-1.21 0.21634 0.30857 —0.58970 
—1.20 0.21944 0.31148 —0.59321 
—1.19 0.22257 0.31439 —0.59669 
—1.18 0.22572 0.31731 —0. 60015 
0.22891 0.32023 —0. 60358 
—1.16 0.23213 0.32315 —0. 60699 
—1.15 0.23538 0.32608 —0.61037 
—1.14 0.23865 0.32902 —0.61373 
—1.13 0.24196 0.33195 —0.61706 
—1.12 0.24529 0.33489 —0. 62037 
—1.11 0.24865 0.33783 —0.62365 
—1.10 0.25205 0.34078 —0. 62690 
—1.09 0.25547 0.34373 —0. 63013 
—1.08 0.25892 0.34667 —0. 63333 
—1.07 0.26240 0.34963 —0.63650 
—1.06 0.26591 0.35258 —0.63965 
—1.05 0.26945 0.35553 —0.64276 
—1.04 0.27302 0.35849 — 0.64585 
—1.03 0.27662 0.36144 —0.61891 
—1.02 0.28025 0.36440 —0.65191 

01 0.28391 0.36736 —0.65494 
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v 
—1.00 0.28760 0.37031 —0.65791 
—0.99 0.29132 0.37327 —0. 66086 
—0.98 0.29507 0.37623 —0.66377 
—0.97 0.29884 0.37918 —0.66665 
—0.96 0.30265 0.38214 —0.66950 
—0.95 0.30649 0.38509 —0.67232 
—0.94 0.31035 0.38805 —0.67512 
—0.93 0.31425 0.39100 —0.67788 
—0.92 0.31817 0.39395 —0.68060 
—0.91 0.32213 0.39690 —0.68330 
—0.90 0.32611 0.39985 —0.68597 
—0.89 0.33012 0.40279 —0.68860 
—0.88 0.33416 0.40573 —0.69121 
—0.87 0.33824 0.40867 —0.69378 
—0.86 0.34234 0.41161 —0.69632 
—0.85 0.34647 0.41454 —0.69883 
—0.84 0.35063 0.41747 —0.70130 
—0.83 0.35482 0.42040 —0.70375 
—0.82 0.35904 0.42332 —0.70616 
—0.81 0.36328 0.42624 —0.70854 
—0.80 0.36756 0.42915 —0.71088 
—0.79 0.37187 0.43206 —0.71320 
—0.78 0.37620 0.43497 —0.71548 
—0.77 0.38057 0.43787 —0.7177 
—0.76 0.38496 0.44076 —0.71994 
—0.75 0.38938 0.44365 —0.72212 
—0.74 0.39383 0.44654 —0.72427 
—0.73 0.39831 0.44942 —0.72639 
—0.72 0.40282 0.45230 —0.72848 
—0.71 0.40736 0.45517 —0.73053 
—0.70 0.41192 0.45803 —0.73255 
—0.69 0.41652 0.46089 —0.73453 
—0.68 0.42114 0.46374 —0.73648 
—0.67 0.42579 0.46658 —0.73840 
—0.66 0.43047 0.46942 —0.74029 
—0.65 0.43518 0.47225 —0.74215 
—0.64 0.43992 0.47508 —0.74397 
—0.63 0.44468 0.47789 —0.74576 
—0.62 0.44948 0.48071 —0.74751 
—0.61 0.45430 0.48351 —-0.74924 
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n v 
—0.60 0.45915 0.48630 —0.75093 
—0.59 0.46402 0.48909 —0.75259 
—0.58 0.46893 0.49187 —0.75422 
—0.57 0.47386 0.49465 —0.75581 
—0.56 0.47882 0.49741 —0.75737 
—0.55 0.48381 0.50017 —0.75890 
—0.54 0.48883 0.50292 —0.76040 
—0.53 0.49387 0.50566 —0.76187 
—0.52 0.49894 0.50839 —0.76330 
—0.51 0.50404 0.51111 —0.76470 
—0.50 0.50916 0.51382 —0.76607 
—0.49 0.51431 0.51653 —0.76741 
—0.48 0.51949 0.51923 —0.76872 
—0.47 0.52470 0.52191 —0.77000 
—0.46 0.52993 0.52459 —0.77124 
—0.45 0.53519 0.52726 —0.77246 
—0.44 0.54047 0.52992 —0.77364 
—0.43 0.54579 0.53257 —0.77479 
—0.42 0.55113 0.53521 —0.77592 
—0.41 0.55649 0.53784 —0.77701 
—0.40 0.56188 0.54047 —0.77807 
—0.39 0.56730 0.54308 —0.77910 
—0.38 0.57274 0.54568 —0.78010 
—0.37 0.57821 0.54827 —0.78107 
—0.36 0.58371 0.55085 —0.78202 
—0.35 0.58923 0.55342 —0.78293 
—0.34 0.59478 0.55599 —0.78381 
—0.33 0.60035 0.55854 —0.78467 
—0.32 0.60595 0.56108 —0.78549 
—0.31 0.61157 0.56361 —0.78629 
—0.30 0.61722 0.56613 —0.78706 
—0.29 0.62289 0.56864 —0.78780 
—0.28 0.62859 0.57114 —0.78851 
—0.27 0.63432 0.57362 —0.78920 
—0.26 0.64007 0.57610 —0.78985 
—0.25 0.64584 0.57857 —0.79018 
—0.24 0.65164 0.58102 —0.79108 
—0.23 0.65746 0.58347 —0.79166 
—0.22 0.66331 0.58590 —0.79221 
0.66918 0.58833 —0.79273 
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—0.20 0.67507 | 0.59074 —0).79322 
—0.19 0.68099 0.59314 —0.79369 
—0.18 0.68694 | 0.59553 —0.79413 
—0.17 0.69290 | 0.59791 —0.79455 
—0.16 0.69889 0.60028 —0.79494 
—0.15 0.70491 0.60263 —0.79530 
—0.14 0.71095 0.60498 —0.79564 
—0.13 0.71701 0.60731 —0.79596 
—0.12 0.72309 0.60963 —0.79625 
—0.11 0.72920 0.61195 —0.79651 
—0.10 0.73533 0.61425 | —0.79676 
—0.09 0.74149 0.61653 | 0.79697 
—0.08 0.74766 0.61881 | —0.79717 
—0.07 0.75386 0.62108 | 0.79734 
—0.06 0.76008 0.62333 —0.79748 
—0.05 0.76633 0.62558 —0.79761 
—0.04 0.77260 0.62781 —0.79771 
—0.03 0.77888 0.63003 —0.79779 
—0.02 0.78520 0.63224 —0.79784 
—0.01 0.79153 0.63443 —0.79787 
0.00 0.79788 0.63662 —0.79788 
0.01 0.80426 0.63879 —0.79787 
0.02 0.81066 0.64096 —0.79784 
0.03 0.81708 0.64311 —0.79779 
0.04 0.82352 0.64525 —0.79771 
0.05 0.82999 0.64738 —0.79762 
0.06 0.83647 0.64949 —0.79750 
0.07 0.84298 0.65160 —0.79736 
0.08 0.84950 0.65369 —0.79721 
0.09 0.85605 0.65578 —0.79703 
0.10 0.86262 0.65785 —0.79683 
0.11 0.86921 0.65991 —0.79662 
0.12 0.87582 0.66196 —0.79638 
0.13 0.88245 0.66399 —0.79613 
0.14 0.88910 0.66602 —0.79585 
0.15 0.89577 0. 66803 —0.79556 
0.16 0.90216 0.67063 —0.79525 
0.17 0.90917 0 6720: —0.79492 
0.18 0.91590 | 0.67400 —0.79458 
0.19 0.92265 0.67597 —0.79421 
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n v 

0.20 | 9.92942 0.67793 —0.79383 
0.21 | 0.93620 0.67988 —0.79343 
0.22 0.94301 0.68181 —0.79301 
0.23 0.94984 0.68373 ~0.79258 
0.24 0.95669 0.68565 —0.79213 
0.25 0.96355 0.68755 —0.79167 
0.26 0.97044 0.68944 0.79119 
0.27 0.97734 0.69132 —0.79069 
0.28 0.98427 0.69318 —0.79017 
0.29 0.99121 0.69504 —0.78964 
0.30 0.99817 0.69689 —0.78910 
0.31 1.00514 0.69872 —0.78854 
0.32 1.01214 0.70054 ~0.78797 
0.33 1.01915 0.70236 —0.78738 
0.34 1.02619 0.70416 —0.78677 
0.35 1.03324 0.70595 —0.78616 
0.36 1.04031 0.70773 —0.78552 
0.37 1.04739 0.70950 —0.78488 
0.38 1.05450 0.71125 —0.78422 
0.39 1.06162 0.71300 —0.78355 
0.40 1.06876 0.71474 —0.78286 
0.41 1.07591 0.71646 —0.78216 
0.42 1.08309 0.71818 —0.78145 
0.43 1.09028 0.71988 —0.78073 
0.44 1.09748 0.72158 —0.77999 
0.45 1.10471 0.72326 —0.77924 
0.46 1.11195 0.72493 —0.77848 
0.47 1.11921 0.72659 —0.77771 
0.48 1.12648 0.72825 —0.77692 
0.49 1.13377 0.72989 —0.77613 
0.50 1.14108 0.73152 —0.77532 
0.51 1.14840 0.73314 —0.77450 
1.15574 0.73475 —0.77367 
0.53 1.16310 0.73635 —0.77283 
0.54 1.17047 0.73794 —0.77198 
0.55 1.17785 0.73952 ~0.77112 
0.56 1.18526 0.74109 —0.77025 
0.57 1.19268 | 0.74265 —0.76937 
0.58 1.20011 0.74420 —0. 76847 
0.59 1.20756 | 0.74574 —0.76757 
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v 
0.60 1.21503 0.74727 —0.76666 
0.61 1.22251 0.74879 —0.76574 
0.62 1.23000 0.75030 —0.76481 
0.63 1.23751 0.75180 —0.76388 
0.64 1.24504 0.75329 —0.76293 
0.65 1.25258 0.75478 —0.76197 
0.66 1.26013 0.75625 —0.76101 
0.67 1.26770 0.75771 —0.76004 
0.68 1.27529 0.75916 —0.75906 
0.69 1.28289 0.76061 —0.75807 
0.70 1.29050 0.76204 —0.75707 
0.71 1.29813 0.76346 —0.75607 
0.72 1.30577 0.76488 —0.75506 
0.73 1.31342 0.76628 —0.75404 
0.74 1.32109 0.76768 —0.75301 
0.75 1.32878 0.76907 —0.75198 
0.76 1.33648 0.77044 —0.75094 
0.77 1.34419 0.77181 —0.74989 
0.78 1.35191 0.77317 —0.74884 
0.79 1.35965 0.77453 —0.74778 
0.80 1.36740 0.77587 —0.74671 
0.81 1.37517 0.77720 |. —0.74564 
0.82 1.38295 0.77852 —0.74456 
0.83 1.39074 0.77984 —0.74347 
0.84 1.39854 0.78115 —0.74238 
0.85 1.40636 0.78244 —0.74128 
0.86 1.41419 0.78373 —0.74018 
0.87 1.42204 0.78501 —0.73907 
0.88 1.42989 0.78629 —0.73796 
0.89 1.43776 0.78755 —0.73684 
0.90 1.44564 0.78881 —0.73572 
0.91 1.45354 0.79005 —0.73459 
0.92 1.46144 0.79129 —0.73346 
0.93 1.46936 0.79252 —0.73232 
0.94 1.47729 0.79374 —0.73118 
0.95 1.48524 0.79496 —0.7 
0.96 1.49319 0.79616 —0.72888 
0.97 . 1.50116 0.79736 —0.72772 
0.98 1.50914 0.79855 —0.72656 
0.99 1.51713 0.79973 —0.72540 
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0.80090 —0.72423 
0.80207 | 0.72306 
0.80323 | 0.72189 
0.80137 —0.72071 
0.80552 —0.71953 
0.80665 —0.71834 
0.80778 —0.71715 
0.80890 —0.71596 
0.81001 —0.71477 
0.81111 —0.71357 
0.81221 —0.71237 
0.81330  —0.71117 
0.81438 | —0.70996 
0.81545 | 0.70875 
0.81652 | 0.70754 
0.81758 —0.70633 
0.81863 —0.70511 
0.81968 —0.70389 
0.82071 —0.70267 
0.82175 —0.70145 
0.82277 —0.70023 
0.82379 —0.69900 
0.82480 —0.69778 
0.82580 —0.69655 
0.82680 —0.69532 
0.82779 —0.69408 
0.82877 —0.69285 
0.82975 —0.69162 
0.83071 —0.69038 
0.83168 —0.68914 
0.83263 —0.68790 
0.83358 —0. 68666 
0.83453 —0.68542 
0.83546 —0.68418 
0.83640 —0.68294 
0.83732 —0.68170 
0.83824 —0.68045 
0.83915 —0.67921 

24 0.84006 —0.67796 
0.84095 —0.67672 


4 
1.00 | 1.5251 
: 1.05 | 1.5653 
1.07 | 1.5814 
1.09 1.5976 
| 
110 | 1.6058 
| 
1.13 1.630 
1.14 | 1. 638: 
1.16 | 1.6547 
1.19 1.679 
1.24 | 1.7208 
| 
1.250 
1.26 | 1.737 
ee 1.27 1.745; 
1.28 1.753 
1.30 | 1.770 
1.31 | 1.778 
1.32 | 1.787 
1.33 1.795 
pea 1.34 1.803 
| 
1.36 1.820 
1.38 1.837 
4 1.39 | 1.84 
| 
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n | | ¢ 

1.40 1.85106 0.84185 | 0.67547 
141 1.86248 0.81274 0.67122 
1.42 1.87091 0.81362 —0.67298 
1.43 1.87935 0.84449 —0.67173 
1.41 1.88780 0.81536 —0.67048 

| 

1.45 1.89626 0.84622  —0.66923 
1.46 1.90473 0.81708 0.66799 
1.47 1.91320 0.84793 | —0.66674 
1.48 1.92168 0.84878 | 0.66549 
1.49 1.93018 0.84962 | —0.66424 
1.50 1.93868 9.85045 0.66300 
1.51 1.94719 0.85128  =0.66175 
1.52 1.95570 0.85211 —0.66050 
1.53 1.96423 0.85292 —0.65926 
1.54 1.97276 0.85373 —0.65801 
1.55 1.98130 0.85454 —0.65676 
1.56 1.98985 0.85534 —0.65552 
1.57 1.99841 0.85614 —0.65427 
1.58 2.00697 0.85693 —0.65303 
1.59 2.01555 0.85771 —0.65179 
1.60 2.02413 0.858 19 | —0.65054 
1.61 2.03272 0.85927  —0.61930 
1.62 2.04131 | 0.86603 —0. 64806 
1.63 2.04992 0.86080 —0.64682 
1.64 2.05853 0.86156 —0.64558 
1.65 2.06715 0.86231 —0.64434 
1.66 2.07578 0.86306 —0.64310 
1.67 2.08441 0.86380 —0.64186 
1.68 2.09305 0.86454 —0.64062 
1.69 2.10170 0.86527 —0.63939 
1.70 2.11036 0.86600 —0.63815 
1.71 2.11902 0.86673 —0. 63692 
1.7 2.12769 0.86744 0.63569 
1.73 2.13637 0.86816 —0. 63446 
1.74 2.14506 0.86887 —0. 63323 
1.75 2.15375 0.86957 —0.63200 
1.76 2.16245 0.87027 —0.63077 
1.77 2.17115 0.87096 —0. 62955 
1.7 2.17987 0.87165 | 62832 
1.79 2.18859 0.87234 | —0.62710 
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n v | 
1.80 2.19731 | 0.87302 —0. 62588 
1.81 2.20605 0.87370 —0. 62466 
1.82 2.21479 | 0.87437 —0.62344 
1.83 2.22353 0.87504 —0.62222 
1.84 2.23229 | 0.87570 —0.62100 
1.85 2.24105 0.87636 —0.61979 
1.86 2.24981 0.87701 —0.61857 
1.87 2.25859 0.87766 —0.61736 
1.88 2.26737 0.87831 —0.61615 
1.89 2.27615 0.87895 —0.61494 
1.90 2.28495 0.87958 —0.61374 
1.91 2.29375 0.88022 —0.61253 
1.92 2.30255 ! 0.88084 —0.61133 
1.93 2.31136 0.88147 —0.61013 
1.94 2.32018 0.88209 —0. 60893 
1.95 2.32900 0.88270 —0.60773 
1.96 2.33783 0.88331 —0. 60654 
1.97 2.34667 0.88392 —0. 60534 
1.98 2.35551 0.88453 —0.60415 
1.99 2.36436 0.88513 —0. 60296 
2.00 2.37322 0.88572 —0.60177 
2.01 2.38208 0.88631 —0. 60059 
2.02 2.39094 0.88690 —0.59940 
2.03 2.39981 0.88748 —0.59822 
2.04 2.40869 0.88806 —0.59704 
2.05 2.41757 0.88864 —0.59586 
2.06 2.42646 0.88921 —0.59469 
2.07 2.43536 0.88978 —0.59351 
2.08 2.44426 0.89035 —0.59234 
2.09 2.45317 0.89091 —0.59117 
2.10 2.46208 0.89146 —0.59000 
2.11 2.47100 0.89202 —0. 58884 
2.12 2.47992 0. 89257 —0.58767 
2.13 2.48885 0.89311 —0.58651 
2.14 2.49778 0.89366 —0.58535 
2.15 2.50672 0.89420 —0.58420 
2.16 2.51566 0.89473 —0.58304 
2.17 2.52461 0.89526 —0.58189 
2.18 2.53357 0.89579 —0.58074 
2.54253 0.89632 —0.57959 
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2.20 2.55150 0.89684 —0.57845 
2.21 2.56047 0.89736 —0.57730 
2:22 2.56944 0.89787 —0.57616 
2.23 2.57842 0.89839 —0.57502 
2.24 2.58741 0.89890 —0.57389 
2.25 2.59640 0.89940 —0.57275 
2.26 2.60540 0.89990 —0.57162 
2.27 2.61440 0.90040 —0.57049 
2.28 2.62341 0.90090 —0.56936 
2.29 2.63242 0.90139 —0.56824 
2.30 2.64143 0.90188 —0.56712 
2.31 2.65046 0.90236 —0.56600 
2.32 2.65948 0.90285 —0.56488 
2.33 2.66851 0.90333 —0.56376 
2.34 2.67755 0.90380 —0.56265 
2.35 2.68659 0.90428 —0.56154 
2.36 2.69563 0.90475 —0.56043 
2.37 2.70468 0.90521 —0.55933 
2.38 2.71374 0.90568 —0.55822 
2.39 2.72280 0.90614 —0.55712 
2.40 2.73186 0.90660 —0.55602 
2.41 2.74093 0.90705 —0.55493 
2.42 2.75000 | 0.90751 —0.55384 
2.43 2.75908 0.90796 —0.55274 
2.44 2.76816 0.90840 —0.55166 
2.45 2.77725 0.90885 —0.55057 
2.46 2.78634 0.90929 —0.54949 
2.47 2.79543 0.90973 —0.54841 
2.48 2.80453 0.91016 —0.54733 
2.49 2.81364 0.91060 —0.54625 
2.50 2.82274 | 0.91103 —0.54518 
2.51 2.83186 0.91145 —0.54411 
2.52 2.84097 0.91188 —0.51304 
2.53 2.85009 0.91230 —0.54197 
2.54 2.85922 0.91272 —0.54091 
2.55 2.86835 0.91314 —0.53985 
2.56 2.87748 0.91355 | 0.53879 
2.57 2.88662 | 0.91396 —0.53774 
2.58 2.89576 | 0.91437 | 0.53668 
2.59 2.90491 0.91478 | 0.53564 
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2.91406 
2.92321 
2.93237 
2.94153 
2.95070 


2.95987 
2.96904 
2.97822 
2.98740 
2.99658 


3.00577 
3.01497 
3.02416 
3.03336 
3.04257 


3.23667 
3.24595 
3.25523 
3.26452 
3.27381 


0.91518 
0.91558 
0.91598 
0.91638 
0.91677 


0.91716 
0.91755 
0.91794 
0.91832 
0.91870 


0.91908 
0.91946 
0.91983 
0.92021 
0.92058 


0.92094 
0.92131 
0.92167 
0.92203 
0.92239 


0.92275 
0.92310 
0.92346 
0.92381 
0.92415 


0.92450 
0.92484 
0.92518 
0.92552 
0.92586 


0.92620 
0.92653 
0.92686 
0.92710 
0.92752 


0.92784 
0.92817 
0.92849 
0.92881 
0.92913 


—0.53459 
—0.53354 
—0.53250 
—0.53146 
—0.53042 


—0.52939 
—0.52835 
—0.52732 
—0.52630 
—0.52527 


—0.52425 
—0.52323 
—0.52221 
—0.52120 
—0.52018 


—0.51918 
—0.51817 
—0.51717 
—0.51616 
—0.51516 


—0.51417 
—0.51318 
—0.51218 
—0.51119 
—0.51021 


—0.50923 
—0.50824 
—0.50726 
—0.50629 
—0 50532 


—0.50434 
—0.50338 
—0.50241 
—0.50145 
—0.50049 


—0.49953 
—0.49857 
—0.49762 
—0 49667 
19572 
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| 
| 2. | | | 
2. 
2. 
2. | 
| 
2.66 
2.67 | 
2.68 
2.69 | 
2.70 
2.71 
2.72 
2.73 
2.74 
2.76 3.06098 
2.78 3.07942 
= 3.09787 
| 3.10710 
3.11638 
| 3.14405 
3.15329 
3.16254 
= 
3.17180 
3.18105 
3.19031 
| 
| 3.19958 
| 3.20885 | | 
| 3.21812 
i | 3.22739 
| 
| | | 
| 
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n v 
3.0 3.28310 0.92044 | —0.49178 
3.37620 0.93250 —0.48546 
3:2 3.46959 —0.47640 
3.3 3.56327 0.93809 | —0.46758 
3.4 3.65720 0.94065 0.45900 

3.5 3.75139 0.94307 | —0.45066 
3.6 3.84581 0.94535 —0. 44255 

3.94046 0.94751 —(0. 43466 
3.8 4.03531 0.91956 —0.42700 
3.9 4.13037 0.95149 —0.41954 
4.0 4.22561 0.95333 41230 
4.1 4.32103 0.95507 —0. 40526 
4.2 4.41662 0.95672 —0.39841 
4.3 4.51237 0.95828 —0.39175 
4.4 4.60827 0.95977 —0.38528 
4.5 4.70432 0.96119 —0.37898 
4.6 4.80051 0.96253 —0.37286 
4.7 4.89682 0.96381 —0.36691 
4.8 4.99327 0.96503 —0.36111 
4.9 5.08983 0.96620 —0.35547 
5.0 5.18650 0.96730 —0.34998 
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STATISTICAL ANALYSIS OF SYMMETRICAL SLOPE-RATIO 
ASSAYS OF ANY NUMBER OF TEST PREPARATIONS 


M. CLARKE 


National Institute for Research in Dairying, 
Shinfield, Berkshire 


1. INTRODUCTION 


ee statistical analysis of slope-ratio assays, in which the response 
is linearly related to a power of the dose, has been developed for 
assays of one test preparation only by Finney (6, 7) & Wood (9). Bliss 
(1) has presented a method of analysis for a slope-ratio assay of any 
number of test preparations, and more recently (2) he has described 
routine methods of computation for such assays. 

An alternative method of analysis for a symmetrical slope-ratio 
assay of any number of test preparations is given here: the computation 
required, particularly in testing validity, is shorter than that described 
by Bliss (2), and although the method is of most use when more than 
one test preparation is assayed, its application slightly simplifies the 
analysis given by Finney (7) for a symmetrical assay of only one test 
preparation. 

The optimum deployment of material and effort in slope-ratio assays 
has been discussed elsewhere (3, 10). Common-zero assays, that is, 
assays with tests at zero dose or a basic dose and at higher doses for 
each preparation, have been recommended, and the arguments in 
favour of a common-zero design are strengthened when the number of 
test preparations to be assayed is increased. The method of analysis 
for symmetrical assays without a common zero dose will be described 
here briefly, but attention will be directed chiefly to assays of the 
common-zero type. 

The method of analysis developed here refers to assays in which the 
response to each preparation is linearly related to the dose itself; its 
adaptation to cases where the response is proportional to some power, 
other than unity, of the dose depends merely on a suitable choice of 
dose-levels. Consideration will be confined only to the cases of assays 
with equal numbers of dose-levels for all preparations, equal spacing 
between consecutive doses of each preparation and equal replication 
for all treatments. For a given number of test subjects some increase 


370 


: 
| 
ag 
— 
| WA 


SLOPE-RATIO ASSAYS 371 


in precision results from having more replication for the standard 
preparation; the increased complexity of the computations, however, 
may offen outweigh the advantages of such an arrangement, and if 
there are to be any block restrictions, equal replication of all treat- 
ments will probably be preferable. The method of analysis described 
by Blizs (1, 2) does allow for higher replication of the zero dose and 
standard preparation than of the test preparation; adaptation of the 
method deseribed in this paper to meet this case introduces some 
simplification, but for brevity is not presented here. 

The form of analysis is first illustrated by means of an example. 
and general formulae are given in later sections. 


2. NOTATION 


In the general treatment of assays Invelving more than two prepara- 
tions, the concise description of the designs may cause a slight difficulty: 
for example, 2 7-point assay might be an assay with a zero dose and 2 
preparations each at 3 levels, or an assay with uw zero dose and 3 prepara- 
tions each at 2 levels. To avoid this type of uncertainty, the following 
notation for symmeirical assays is proposed. An assay with a zero 
dose and v preparations each at /& levels will be referred to as a {1 + o(k)}- 
point assay, and an assay with no common zero dose and v preparations 
each at k levels as a {v(i')}-point assay. In this notation the first of 
the cases mentioned above would he described as a {1 + 2(3)}-point 
assay, and the second as a {1 + 3(2)}-point assay. 


3. METHOD OF ANALYSIS 
Example: a \1 + 3(2)}-point assay 


The data used for an example are taken from the results of an assay 
of vitamin 2, in the dried contents of rumina of calves. Although 9 
preparations were in fact tested, each at 3 dose levels, for present 
purposes attention has been confined to 3 preparations each at 2 levels. 
Since there were only 2 replications, restriction of the analysis to re- 
sponses to 7 treatments gives only 7 degrees of freedom for the estimation 
of error variation, but the example will serve to show the principles of 
the analysis. 

Table 1 below shows the individual observations y, quoted in units 
of 10° X optical density. The dose scales for each preparation have 
been standardized so that the higher dose of each preparation corresponds 
to Dunit; with equal spacing of the doses, the lower dose of each prepara- 
tion then corresponds to 4 unit. Phe method of computation and use 
of the quantities Q, band 7/7 will be explained. 
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TABLE 1.) ASSAY OF VITAMIN Be IN DRIED RUMEN CONTENTS: RESPONSES AND 
COMPUTATION TABLE. (FIGURES DUE TO J. E. FORD) 

1 unit of each preparation is equal to 1 ml. 

1 unit of standard preparation contains 0.01 we. pyridoxine 

Dilution of test preparation (1): 1,490 g. ‘ml. 

Dilution of test preparation (2): 1493 ¢. ‘ml. 

Responses are given in 103 X optical density: treatment totals shown in brackets 


Dose (units) | Control Standard Test 1 Test 2 Total 
0 95 95 | 
(190) 
190 190 190 185 203 200 
(880) (375) (403) 
1 267 270 270 262 280 268 
| (537) | (532) (548) 
| | 
Total ' 190 O17 907 951 2965 
Q | 1283 1178 1598 4059 
b i 167.125 164.125 176.125 
H 223 218 258 699 


For each preparation together with the ‘blank’ or zero dose, the 
mean responses are approximately linearly related to the dose. Straight 
lines, intersecting at zero dose, are fitted to the mean responses, by 
estimating the multiple linear regression equation 


where Y is the estimated mean response to 2s , X71 , ro , the doses of 
the standard and test preparation respectively, and bs , br, , bro are 
estimates of the mean increases in response per unit increase in the 
dose of the standard and test preparations respectively. 

The estimated potencies of the test preparations relative to the 
standard are then b7,/bs and bDys/bs. 

The symmetry of the design permits the regression coefficients and 
the sum of squares accounted for by the regression to be obtained by 
direct substitution in simple formulae. This avoids the necessity of 
solving a set of simultaneous linear equations, and therefore has ad- 
vantages, which obtain even with assavs involving only two preparations. 
Incidental computational advantages appear with this method, in that 
it is not generally necessary to carry a large number of figures in inter- 
mediate computations to avoid errors due to rounding-olf. 

A value Q is first computed for each preparation, where, for example. 
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Here 3S, and WS, denote the total responses to the lower and higher 
doses respectively of the standard preparation, and G denotes the 
total of the responses to all treatments. “The values of the regression 
coefficients are then immediate as, for example, 


where nis the number of replications for cach treatment, 


Uere, 


Qs = 7(380 + 1074) — S895 = 1288 


S(Q) = 4059 and = $67,125. 


In order to test the significance of departures from the mathematical 
model and from the underlying biological assumptions, a quantity H 
is calculated for cach preparation, giving the total response to zero 
dose as estimated by a straight regression line fitted to the mean re- 
sponses to non-zero doses of that preparation. These quantities H are 
calculated as 


H = 28, 


and similarly for cach test preparation. 

The assumptions made in this type of biological assay have been 
discussed by Jerne & Wood (8). If the assay is statistically valid, the 
results will be consistent with the hypothesis of a linear relationship 
between response and dose throughout the range of responses observed. 
Accordingly the difference between each value of JZ and that of 7, , the 
total response to zero dose, may be tested for significance; it is generally, 
however, more convenient to test the mean value, over all preparations, 
of this difference. 

The fundamental validity (8) of the assay requires that the values 
of IT for the different preparations should, apart from ehance variation, 
he equal, since regression lines fitted separately to the responses for 
each preparation should intersect at zero dose. 

These tests are conveniently incorporated into an analysis of vari- 
ance, and the results for this example are summarized in Table 2 
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TABLE 2. ANALYSIS OF VARIANCE FOR THE {1 4 3(2)}-POINT ASSAY OF VITAMIN I, 


Sourees of Degrees of Sums of Mean 
variation freedom squares squares 
Correction for mean 627 ,944.64 
Between treatments 6 49 670.86 8,278.48 
Regression 3 49,229.17 16,409.72 
Blanks 1 346.69 346.69 
Intersections 2 95.00 47.50 
Within treatments 7 125.50 17.93 
Total 13 49 ,796 .36 


Sums of squares for differences hetween treatment means, with 6 degrees 
of freedom, and for differences between responses to the same treatment, 
with 7(n — 1) degrees of freedom, are calculated in the usual way. 
That part of the sum of squares accounted for by the regression, with 
3 degrees of freedom, is easily obtained as 


(bsQs + briQri + br2Qr2)/14. 


In this example the sum of squares for regression is therefore [(167.125 X 
1283) + (164.125 X 1178) + (176.125 & 1598)]/14, ie. 49229.17. Since 
the values of b are used as multipliers in this way, it is advisable to 
calculate them to about 6 significant figures. 

To test the validity of the linear regression model a component, 
with 1 degree of freedom, corresponding to Finncy’s term Blanks (7), 
is calculated as 


[S(H) — 37,]?/24n = (699 — 570)°/48 = 346.69 


For a test of the fundamental validity of the assay, the component 
of the sum of squares due to differences between the values of JJ for 
different preparations, is calculated as 


SUH") /5n — /1in 
= |(223)” + (218)* + (258)*]/10 -- (699)7/30 = 95.00. 


This sum of squares has 2 degrees of freedom, and corresponds to 
Finney’s Intersections (7). A useful check of the computations is ob- 
tained here since the total of the sums of squares for regression, blanks 
and intersections should exactly equal the sum of squares for differences 
between treatment means. 

Comparison of the mean square for intersections with the error mean 
square shows that the observed responses are not inconsistent with the 
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hypothesis that straighé regression Imes through the responses to 
non-zero doses intersect at zero dose. ‘The mean square for blanks, 
however, is highly significant, and the validity of the statistical model 
is therefore questionable. ‘There is not, on the basis of these results 
alone, sufficient evidence to examine the adequacy of a linear regression 
model over a restricted part of the observed range, such as for the 
responses to non-zero doses only, and the results of this example will 
not be used to obtain estimates of relative potencies. 

If the evidence supports the foundations of the assay, estimates 
of the relative potencies and of their fiducial limits may usefully be 
made. The relative potencies of the test prepurations are given by 


R, = br,/ds and 


so that one unit of the first test preparation is estimated to be as effective 
as Ft, units of the standard preparation, and similarly for the other 
test preparation. These values are easily converted into terms of the 
original units. 

The estimation of the fiducial limits of the relative potencies pre- 
sents little difficulty. If bs is sufficiently larger than its standard error, 
approximate estimates of the fiducial limits may be used; in most 
microbiological assays this condition will be satisfied, but the point 
may easily be tested by computing the value of g as 


g = 


where s is the error mean square and the value of ¢ appropriate to 
the level of probability and number of degrees of freedom for error is 
found from standard tables. If then g is found to be of the order of 
0.1 or less (7), the fiducial limits of any one relative potency R may 
be calculated as R + tV V(R), where the variance of R is given by 


V(R) = s*[17(R — 1)? + 16R]/10nb? 


If the value of g is not sufficiently small, the exact fiducial limits are 
obtained from Fieller’s formula (5) as 


[ ral + 2R 17 - (1 q) 


The {1 + v(k)}-point assay 


The method of analysis described above may easily be extended for 
symmetrical assays with v preparations cach at / levels. 

The only new feature of the analysis when / is greater than 2 is that 
when the sums of squares for regression, blanks and intersections have 
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been subtracted from the sum of squares for differences between treat- 
ments, a component with v(f — 2) degrees of freedom is left, due to 
deviations from straight regression lines fitted to the responses to non- 
zero doses. This component may, if desired, be subdivided into quad- 
ratic, cubic, etc. components for each preparation separately, using the 
method of orthogonal coefficients (see, for example, 3 & 4). 

The formalae required in the analysis of a {1 + v(k)}-point assay 
ure summarized in ‘Table 3; their use will be clear from the preceding 
section. These formulae are based on equally spaced dose-levels 
QO, 2,k, , 1 for each preparation. The estimated values of the 
slopes of the dose-response lines depend on the scale used for the dose- 
levels: the formulae presented here will, for example, differ by a factor 
i: from those given by Bliss (2), which are based on coded dose-levels 
0, 1,2, --- ,k. The difference does not, of course, effect the estimate 
of relative potency, but the formulae given by Bliss for the estimated 
variance and fiducial limits of the relative potency are not applicable 
to the dose-scale used in this paper, and conversely. The scale chosen 
in this paper has slight computational advantages, and is the same as 
that used by Finney in his 1951 paper (7). 

The formula for H im any particular design may be written down 
easily by noting that the coefficients of the treatment totals are suc- 
cessively reduced by 3, with the check that the coefficient of the total 
for the highest dose is (—3) times the coefficient of the total for the 
lowest dose. In general H is not equal to the total response to zero 
dose as estimated by a straight line fitted to the responses to non-zero 
doses, but is equal to k(k — 1)/2 times that value. The exact fiducial 
limits of R are given in Table 3; in most microbiological assays, however, 
g is found to be negligible, so that considerable simplification is possible, 
as also shown in Table 3, and with experience of an assay it is usually 
unnecessary to calculate g. 


The {v(k)}-point assay 


It may sometimes be convenient to use designs without a common 
zero dose, and the method of analysis for such assays is similar to that 
already described. The formulae required for making validity tests 
and estimates of relative potencies and their fiducial limits are given 
in Table 3. At least 3 dose-levels for cach preparation are necessary 
for tests of the significance of deviations from linearity of the dose- 
response curve to be possible. 

The formulae in Table 3 apply to standardized dose scales so defined 
that the highest dose of each preparation corresponds to | unit and the 
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lowest. dose to Lok unit, with the intervening doses equally spaced. 
The analysis of variance differs slightly from that for the common-zero 
type of assay in having no component for blanks. 


4. SUMMARY 


A method of analysis for symmetrical slope-ratio assays of any 
number of preparations is presented. It includes specific tests of the 
validity of the assumptions on which the assays are based, and is simpler 
than previously published methods. Some advantage is gained when 
the method is applied to cases where one test preparation only is assayed. 
The calculations required are illustrated with an example, and general 
methods of analysis are described for symmetrical assays with and 
without a common zero dose. 

I should like to thank Mr. C. P. Cox for his encouragement and 
criticism, and Mr. J. E. Ford for permission to use his figures in the 
example. 
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AUTOMATIC MACHINE METHOD OF CALCULATING 
CONTINGENCY X° 


JOHN Skory 


Stalisticul Research Section, American Cancer Society, New York 


N THE September 1951 issue of Biometrics, Leslie (1) presented a hand 
I machine method of calculating X* for an R X C contingency table. 
The method is essentially a reciprocal-multiplication technique. The 
direct or the division method for an automatic machine* is presented 
and compared in this paper. 

In comparing the two methods, the 3 X 4 table given by Leslie 
was the subject for the calculations. Given the original data arranged 
in 3 rows and 4 columns, the following operations were performed and 
the quotients cumulated without transferring to paper any intermediate 
values from the machine: 


69 27 _ 99 5766 
38 4+ 184 183 4 91 80.5766 
7 38? 32? 

38 + + 193 + = 19-6230 
70°. 82”. 55? 

38 + t+ 193 + 97. 97-9048 


* The routine as set forth can be performed on a fully automatic machine such as the Friden, 
Marchant, Monroe and perhaps others. 
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These three cumulated quotients were entered in column 6. 


1.050266 


A 
| 76 69 27: «196 80.5766 
S| 7 70 82 58 | 214 | 97.9048 


38 «184183 91; 496 


These sums of quotients were then divided by the respective row 
totals given in column 5 and the resulting quotients cumulated in one 
operation: 


80.5706 


15.6230 97.9048 
196 


———— = 1.050266 


+ O14 


Then: 
X” = (1.050266 — 1)496 = 24.93 


The unique feature of this method is that, given the original table, 
the only numbers that need to be transferred from the machine to 
the work sheet are those in column 6. 

Results of the Time Tests 

The methods were compared in 6 replicates as to length of time 
required to calculate the NX? value. The 3X 4 table was modified 
slightly by adding one to each cell and mixing up the cells cach time 
the table was used. This minimized the carry-over of values by memory 
but still provided a table of approximately equal difficulty for each 
method in each replicate. 

Before the actual testing both methods were practised for some 
time until the technique was familiar to the operator. The methods 
were randomized within each replicate. 

The mean number of minutes for each method is given in table 1 
and the analysis of variance in table 2. 
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TABLE 1 
EXPERIMENTAL DATA 


Replicate in time | | Mean 
Method | 1 2 3 4 5 6 | (minutes) 
Multiplication | 7.9 7.7 6.1 6.7 64 63] 41.1 | 6.85 
Division ; 6.0 4.5 4.4 4.5 4.4 4.0 26.8 LAT 
| 12.9 12.2 10.5 11.2 10.8 10.3 
TABLE 2 
ANALYSIS OF VARIANCE 
ms. F 
Replicates 5 
Regression 1 1.945 1.945 12.2* 
Remainder 4 .690 
Methods 1 17.042 17.042 107 .2** 
Methods X Replicates 5 .793 159 


*significant at the 5% level 
** significant at the 1% level 


Although there were no apparent differences among replicates, 
(5 d.f.), the regression of methods on time was significant at the 5 
percent level because replicates were consecutive periods of approxi- 
mately 30 minutes each and there was an increase in efficiency in the 
calculation of X? as the method became more familiar to the operator. 
A highly significant difference was found between the two methods 
of calculation tested. The direct method took about 2/3 the time for 
calculation required by the reciprocal-multiplication method (Leslie) 
when an automatic machine was employed. 
The greater efficiency is ascribed to: 
(1) less copying of numbers on paper 
(2) less transfer of large numbers from paper to machine and 
vice-versa 
(3) fewer number of major steps in the calculation 
(4) less worry in determining the correct number of decimals 


REFERENCE 


1. Leslie, P. H.: The calculation of X? for an R X C contingency table. 
Biometrics, 7, 283-286, 1951. 
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QUERIES 


Grorce W. Sxepecor, Mditor 


% QUERY: We have come up with a “missing-plot”’ problem which 
we can’t answer satisfactorily. 

In order to find more suitable lures or baits for fruit fly 
control, we placed 12. wind-rotated cireular suspension platforms in 
infested areas, and placed 4 bait traps at random on each platform. 
Rach of these 4 traps carries a different bait. Hence we have 12 in- 
dependent replicates each with the same 4 baits: A, B,C, and D. 

We have recorded the following fruit fly catches by the 4 baits (over 
3-day period): 


Baits 
Reps A B e D | Totals 
1 8 75 41 47 171 
2 5 52 30 | 29 116 
3 58 |; | 129 
4 3 19 2 | 10 52 
5 17 4] 35 | 79 172 
6 15 34 28 13 | 90 
7 44 | 24 12 94 
8 19 56 23 18 116 
10 20 33 263 | 17 96 
ll 7 25 15 12 59 
12 21 66 39 21 147 
Totals | 132 524 325 | 306 1287 


There was a missing value in Rep. No. & for bait “A. When this 
value is calculated by the usual formula, we get an impossible value, 
i.e., less than zero, e.g. 


+ _ (1382) + (12945) _ 
(3)(11) 


How do vou explain it? Have we made an error somewhere? (Some 
years ago we also found a situation in which the use of a missing-plot 
formula gave us an wnreasonably high vield for a missing plot in a field 
test.) 
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You have made no error in the calculation. The value 
ANSWER: you got will give you the correct estimate of the catch 

for bait A and will lead to the correct analysis of variance. 
For the correction of a small bias in the treatment mean square, sec 
Yates’ article in “The Empire Journal of Experimental Agriculture”, 
1:129 (1933) or *Missing-Plot Techniques” by R. L. Anderson, This 
Journal 2:41 (June, 1946). 

The value of X is not intended as an estimate of the missing datum. 
It is merely a number to be put in the empty space so that you can 
conveniently extract the remaining information from the experiment. 

If vou should apply the appropriate least squares technique, you 
would arrive at the same estimate of the bait mean, together with the 
same analysis of variance, without making any assumption about the 
value of the missing datum. The missing plot formula provides a 
number to be inserted in the table so that the ordinary procedures can 
be used, procedures which are shorter as well as more familiar than the 
least squares solution. An interesting and informative statement of the 
problem may be found in Kempthorne’s “The Design and Analysis of 
Experiments’, sections 6.5.3-5, 9.5-6. 

Incidentally, the catches may well be distributed in Poisson dis- 
tributions rather than normal. Also, some of your numbers of flies are 
small. If there are critical decisions to be made, you might find it 
worthwhile to use the square root transformation. But I suspect you 
will reach the same conclusions by either method. 


QUERY: In a5 X 5 latin square experiment for investigating 
97 the effect of spacing on the yield of millet, the pertinent part of 
the analysis of variance is as follows: 


Source of Variation Degrees of Freedom Mean Squire 
Spacings: | 
Linear Regression | 1 3,960 
Deviations 3 65.7 
Experimental Error | 12 1,056 


For regression, F = 3,960/1,056 = 3.70; Fo; = 4.75 


From the above table it is clear that regression is not significant. 
But if one performs the ordinary (-test for testing the regression co- 
efficient, one gets significance. The calculations are given below: 


b = SXY/SX° = —178/40 = —4.45 


| 

bg 
| 


QUERIES 385 


The standard error of b is 1/65.7/(40)(5) = 0.572 
t= —4.45/0.572 = —7.78 
Since (9, = 3.182, the regression is significant. 


The same result is obtained if the sum of squares for regression was 
tested against the Deviation sum of squares, instead of Error sum of 
squares. 


I would like to know which one of these is the correct procedure. 
The correct: procedure depends on what you consider to 
ANSWER: he the estimate of the real error involved. From the 
information given, I think the real error is estimated by 
the experimental error in the latin square experiment. My opinion is 
that the first procedure is correct and that regression is not significant. 

If the first. two lines in the table were the only information to be 
had, then I would use the second procedure. But from the additional 
information in the third line, it seems clear that I should have been 
misled by what appears to be an unusual underestimate of the true 
variance for deviations. If regression is really linear, and if the other 
theoretical conditions are satisfied, then deviations and experimental 
error estimate the same o°. 

Another way of saying all this is that the experiment seems to have 
been designed for testing all effects against experimental error. The 
4 degrees of freedom for spacings were planned to be separated into 
orthogonal portions suitable to the regression turned up, each portion 
to be so tested. It happened that the spacing yields were closer to 
the linear regression than would be expected from the size of the ex- 
perimental error. This should not effect the planned procedure for 
testing. 


ABSTRACTS 
ENAR (Firru AnNvAL MEETING) 


A joint symposium with the American Society of Horticultural Science 
included five papers on plot technique and two on taste testing and quality 
evaluation. These are to be published as a group and in full at an early 
date. 


WALTER C. JACOB AND ARNOLD KLUTE. (Cornell 
199 University). Sampling Soils for Physical and Chemical Proper- 
ties. 


Samples were taken from an organic amendment experiment to 
determine the relative magnitudes of sources of variance for various 
chemical and physical properties of the soil. Purpose of the paper is 
to put on record these variances in the hope that other people would 
publish similar information. The use of such data for developing a 
sampling plan is described. 


REX LEROY HURST. (Cornell University). Effects of Certain 
200 Soil and Climatic Factors on the Yield and Composition of Corn 
—Variability Study. 


Various soil and climatic factors were measured on each of 160 
field plots, arranged in sixteen replications of ten treatments, in an 
attempt first, to determine the sources of variation in field trials, and 
second, to find out if these sources of variation could be eliminated by 
covariance. 

The factors which were measured are as follows: soil moisture and 
its distribution, soil classification, moisture equivalent, fifteen atmos- 
pheres of tension, bulk density, stand, organic matter, total nitrogen, 
pH, nitrates, phosphates, and potassium. 

The effect of soil moisture and its distribution within a season was 
studied using a modification of the Fisher rainfall technique. 

Soil classification was the most important single variable and 
accounted for approximately fifty percent of the within treatment 
variation. Soil moisture and its distribution accounted for approxi- 
mately the same amount of variation and the two sources of variation 
were largely overlapping. 
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Complete covariance analyses, using the most important factors, 
resulted in only small gains when the block size was small. Appreciable 
gains, up to fifty-eight percent, were realized when larger blocks were 
used. Covariance might be used in place of incomplete block designs, 
when the sources of variation are known. 


WALTER D. FOSTER. (Department of Agricultural Bio- 

201 chemistry, University of West Virginia). The Use of the Dis- 
criminant Function in the Assessment of Human Nutritional 
Status. 


Iluman nutritional status often is estimated from three types of 
measurements: (1) physical signs, customarily classified into four 
degrees by the observing physician; (2) blood levels of nutrients as 
determined by the biochemist; and (3) levels of dietary intake as 
estimated by the nutritionist. 

To the extent that the levels of blood and dietary nutrients are 
correlated with the incidence of physical signs are there possible indi- 
cations of simple cause and effect relationships. Correlations among 
these three types of measurements constitute one step towards the 
assessment of human nutritutional status. 

The discriminant function, with a few modifications, offers a statisti- 
cal procedure for examining these correlations and analyzing the relative 
contribution of various nutrients towards the absence or presence of a 
physical sign. It provides estimates for the degree to which a com- 
hination of blood and dietary nutrients explain the incidence of a 
clinical sign, estimates of the magnitude and direction of the association 
of each nutrient (when considered in combination) with a sign and 
offers tests of significance for these estimates. 


W. T. FEDERER. (Cornell University). Additional Treat- 
202 ; 
ments in Factorial Experiments. 


If more information is desired on one or more combinations of factors 
in a factorial experiment or if some combinations have different error 
variances than others additional plots of the combination are included 
in each replicate of the design. The analysis of such experiments, the 
orthogonality of the contrasts, and the variance of contrasts are dis- 
cussed. A “proportional” factorial scheme is proposed for certain 
experimental situations. Also, biased estimates of contrasts are pro- 
posed which are orthogonal and whieh have smaller variance than the 
unbiased estimates. Tn conclusion, it is recommended that the addi- 


*Biometrician for the cooperative project oi Haran Nutritional Status in the Northeasi Region, 
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tional units or plots be utilized by having an additional replicate on all 
combinations unless it is necessary to obtain a separate estimate of the 
variance for the particular combination or combinations. 


203 EARL L. GREEN. (Ohio State University). The Sampling 
Variance of the Mean of a Two-Threshold Quantitative Trait. 


The problem is to find the sampling variance of the mean of a sample 
from a (supposedly) normal population when the sample yields only 
three proportions: the proportion below one threshold, the proportion 
above the first threshold but below a second threshold, and the pro- 
portion above the second threshold. The variance of the mean will be 
useful in comparing quantitative genetic models with the results of 
breeding experiments when the trait in question presents a graded 
trichotomy. 

The mean and variance of the sample may be computed by (1) trans- 
forming the proportions below the first and below the second thresholds 
into normal deviates, (2) rescaling with the distance between the 
thresholds as a “unit’’, and (3) choosing an arbitrary point, such as 
Threshold 1, as origin. If ¢, and ¢, are two normal deviates corre- 
sponding to the proportions p, below threshold 1 and p, below thres- 
hold 2 (p, including p,), and if ¢, is chosen as a new origin, the mean 
% = —t,/(t, — t,) and the standard deviation s, = 1/(t, — t,). The 
sampling variance of the mean may be found by propagating the error 
in p, and p, into ¢, and ¢, , and thence into #. For the case with ¢, as 
origin, 


— 4) - 


where z, and 2, are the ordinates of the normal curve at /, and ¢, and n 
is the sample size. 


MARIANNE E. BERNSTEIN. (Radcliffe College and Biology 
204 Laboratories, Harvard University). Evidence of Genetic 
Variation of the Human Sex Ratio. 


L. M. Terman (1936) devised psychological tests to test the degree 
of masculinity-femininity and found strong correlation between pro- 
fessions and scores obtained on this test. Working with the upper 
social classes in order to minimize the effect of abortions the author 
found strong correlation between the degree of masculinity respective 
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to femininity of either parent as measured by occupation and the sex 
ratio of the offspring. 

In particular men or women engaged in such “masculine” occupa- 
tions as law, politics, business administration, etc. had significantly 
more male offsprings than men or women engaged in such “feminine” 
occupations as art, music, home economics, social work, etc. A genetic 
model is presented to explain the possible action of the genes involved. 


THEODORE W. HORNER. (Department of Experimental 
205 Statistics, North Carolina State College). Genetic Variance 
Expectations When Gene Interactions are Present. 


Knowledge of the expectations of genetic variances and covariances 
when non-allelic gene interaction is present is useful in determining the 
consequences of non-allelic gene interaction (epistasis). A procedure 
is described for deriving these expectations when non-allelic gene inter- 
action is present. This procedure employs the multinomial distribution. 
Attention is given to the mathematical representation (1) a comple- 
mentary model, (2) a duplicate factor model, (3) a multiplicative model, 
and (4) the model proposed by Wright in 1935. The genetic variance 
among means of F; progenies when complementary gene action is pres- 
ent is derived as an example. Some of the results accruing from the 
application of these techniques are listed. 


H. D. LANDAHL, GEORGE KARREMAN. (University of 
206 Chicago). A Physico-Mathematical Theory of Threshold 
Phenomena: Single and Repetitive Discharges in Excitation. 


A simple physico-chemical model previously studied, which could 
account for a number of threshold phenomena, is generalized to account 
for a single discharge with return to the initial state, as well as for 
repetitive discharges. In the model used the type of response is sensitive 
to Ca** ion concentration. For example, for high Ca** concentrations, 
other quantities being unchanged, there is no threshold. For a lower 
concentration, a threshold appears. This is the change necessary to 
pass between a charged and discharged state. At a slightly lower con- 
centration there is a threshold required to produce what corresponds 
to a single discharge with return to the initial conditions. At still lower 
Ca** concentrations there is an undamped oscillation with rather sharp 
potential changes similar to that in repetitive discharges. However, 
if the Ca** concentration is too low, these subside, decreasing in ampli- 
tude and frequency so that the oscillation completely disappears. 
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GEORGE KARREMAN. (University of Chicago). A Problem 
207 in the Mathematical Biology of Blood Circulation: On the 
Reflection of Pressure Waves in the Arterial System. 


The reflection of pressure waves in a fluid enclosed in an elastic 
tube is treated on the basis of simple physical principles. Reflection 
and transmission coefficients for the case of a sudden change in the 
diameter of the tube, in the elastic properties, or in the thickness of 
the wall, have been derived in terms of the change in the radii and in 
the velocities of propagation of pressure waves in the corresponding 
parts of the tube. The importance of the derived formula for the 
reflection coefficient is illustrated by determining the way in which 
the pressure, resulting from the superposition of an advancing sinusoidal 
pressure wave and its reflected wave, depends on the ratio of the radii, 
of the elasticity moduli, and of the thickness of the walls. This is done 
for the case in which these characteristics suddenly change. The 
clinical implications are pointed out. 


JOSEPH ZUBIN, CARNEY LANDIS. (Columbia University). 
208 Group vs. Individual Statistics in the Evaluation of the Effects 
of Chemical Agents on Humans. 


In determining the effect of certain chemical agents (phenobarbital, 
dexedrine, antihistamines, etc.) on tests of efficiency (flicker fusion, 
reaction time, tapping and subjective feeling-tone) an‘analysis of vari- 
ance of the group data gave non-significant results for the various 
agents but an individual analysis of the data for each subject indicated 
that for some subjects the dosage of each agent was above threshold 
while for others it was not. The relationship between the individual 
analysis, the group analysis and the contrast between performance and 
subjective feeling-tone will be discussed. 


209 C. I. BLISS. (Connecticut Agricultural Experiment Station). 
Fitting the Negative Binomial Distribution to Biological Data. 


The value of the negative binomial distribution in analyzing bio- 
logical counts is enhanced by a practicable maximum likelihood estimate 
of the coefficient k, for which I am indebted to Sir Ronald Fisher. Its 
calculation is illustrated numerically and compared with other methods 
of estimation. The error variances of the negative binomial and tests 
for its agreement with a set of observations are described. A review 
of underlying models, with illustrations from both plants and animals, 
suggests that it has the widest applicability of the known distributions 
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for over-dispersion. With the new technique, a single coefficient & can 
be fitted to a series of similar distributions and tested for homogeneity. 


210 CLIFFORD MALONEY. (Camp Detrick, Frederick, 
Maryland). Some Statistical Aspects of Bacterial Plate Counting. 


It is usual in bacterial plate counting to pour duplicate or triplicate 
plates at each of a series of 10-fold dilutions, and to count only the two 
or three plates showing a number of colonies between certain limits— 
say 30 and 300. This paper gives statistical reasons showing that pour- 
ing only one plate at a closer dilution ratio will yield more accurate results 
if the colonies on all countable plates are counted and properly combined. 


DOUGLAS ROBSON. (Cornell University). Interval 
Estimation of Fish Populations by Means of Marked Members. 


A population consisting of a known number m of marked fish and an 
unknown number of wu of unmarked fish is sampled for the purpose of 
estimating u. The number X of marked fish in a sample of size n is 
assumed to follow the hypergeometric distribution and hence is a 
sufficient statistic for estimating u. 90% and 95% confidence interval 
estimates are constructed by forming the critical seis of size .10 and 
.05 in the z-space which zontair. the maximum number of points for each 
possible value of w. The tabulated intery:|ls for m = 10, n = 1(1)40. 
0 <u < 75; m = 15, = 1(1)15, n = 18,0 < u < 75; m = 20, 
n = 1(1)15, n = 20,0 < u < 75 are compared to the 95% “equal- 
tailed” and the 95% “one-tailed” interval estimates by graphing the 
respective power curves. The effect of modifications in the experiment, 
such as a change of sample size or a change of the proportion of marked 
fish in the population, are also evaluated by a comparison of power 
curves. 

Numerical examples of the following binomiai approximation suggest 
that results very close to the hypergeometric “equal-tailed” interval can 
be obtained very simply: tables of confidence intervals for the parara- 
eter p of the binomial distribution are entered with sample size m and 
with x successes; the uppei *ad lower limits, p and p, are then inserted 
into the inequality. 
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Sustaining Members. The Constitution of the Society was amended 
in 1951 to provide for Sustaining Members and implemented by Article 
No. 11 in the Council By-Laws. Under this By-Law a Committee on 
Sustaining Members was named consisting of Edgar B. Carter, Arthur 
C. DeGraff, Melvin W. Green, W. E. Hambourger, Bernard L. Oser, 
and Lloyd C. Miller, Chairman. As a result of their activities, we are 
happy to report that the Society now has the following eight sustaining 
members: Abbott Laboratories; Armour and Company—Research 
Division; Bristol Laboratories Inc.; Collett-Week-Nibecker, Inc.; Eli 
Lilly and Company; Sharp and Dohme, Incorporated; Smith, Kline 
and French Laboratories; and Wyeth Institute of Applied Biochemistry. 
We greatly appreciate their support, which will make it possible for us 
to issue a new directory of the Society early in 1953. 


Belgium and Belgian Congo. Recently the Society named Dr. Leo- 
pold Martin of Brussels National Secretary for Belgium. He has been 
outstandingly successful in bringing the scientists in his country who are 
interested in biometry into the Society. Primarily as a result of his 
activity, the Society currently has in Belgium and the Belgian Congo 43 
members, who are now developing a regional organization. Of the 43 
members, 16 are in academic institutions. Subject matter interests are 
as follows: agronomy 13, forestry and fishery biology 9, horticulture 3, 
medicine 10, veterinary medicine 2, pharmacology 2, and one each in 
bacteriology, chemistry, mathematics and textile engineering. 

Much of this interest may be traced to the attention which biometry 
is receiving in Belgium. Thus Brussels University has created an Insti- 
tute of Statistics and in its School of Medicine a statistical section under 
Professor Olbrechts and Dr. Martin. Instruction in biometry includes 
lectures on cardiovascular biometry by Dr. Martin in the course on 
Circulatory Dynamics. The Center of Agricultural Physical Chemistry 
at Ghent, part of the I.R.S.I.A., has formed a biometric section. The 
Institut Agronomique de |’Etat 4 Gembloux has organized a_post- 
graduate course in biometry. We may look forward to numerous con- 
tributions from the new Société Quetelet, as it is proposed to call the 
Belgian Region. 
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Denmark. Workers in different branches of biology were invited by 
the Danish members of the Society to discuss the application of bio- 
mathematical and statistical methods to problems of experimental 
design and the evaluation of data. Seventy-nine scientists expressed 
an interest and of these nearly 60 attended a meeting at the State 
Serum Institute in Copenhagen in May. The discussions were started 
by N. K. Jerne on problems concerning the determination of the toxicity 
of bacterial toxins. L. Hammarlund presented measurements of the 
fungicidal effect of certain drugs on tomato plants. The discussion on 
these two subjects was introduced by N. W. Bentzon. F. Quaade 
presented extensive data comparing anthropometric data for fat children 
and other children; G. Seidelin introduced the discussion. An attempt 
to determine the age of sea-swallows based on measurements of eggs 
found in nests was presented by F. Ségaard Andersen, and discussed 
by Arne Jensen. Finally a discussion, introduced by. P. S. Ostergaard, 
took place on the comparison of plant hormones in field experiments, 
presented by Sigurd Andersen. G. Rasch summed up the discussions. 
The meeting was very successful and all participants welcomed the 
idea of arranging similar meetings during the autumn. 


ENAR. The fifth annual meeting of ENAR was held at Cornell 
University, Ithaca, N. Y. on September 8-10, 1952, jointly with the 
American Institute of Biological Sciences. At the Regional business 
meeting on September 10, the following officers were named for 1953; 
Vice-President—S. L. Crump; Sceretary-Treasurer—W . T. Federer; mem- 
bers of the Regional Committee, 1953-1955—T. A. Bancroft and H. W. 
Norton. The scientific program consisted of four sessions and a joint 
symposium with the American Society of Horticultural Science (ASHS). 
On Monday, September 8, sessions were held on ‘‘Contributed Papers’ 
under the chairmanship of R. E. Comstock with papers by W. C. Jacob 
and Arnold Klute, R. L. Hurst, W. D. Foster, and W. T. Federer, and 
on “Statistics and Genetics’ under the chairmanship of C. R. Henderson, 
with papers by E. L. Green, Marianne E. Bernstein, T. W. Horner, and 
C. C. Lowe. The session on Tuesday morning, September 9, was en- 
titled “Mathematical Biology” with E. L. Green presiding and papers by 
H. D. Landahl and George Karrerman, George Karreman, and Joseph 
Zubin and Carney Landis. The Tuesday afternoon session was a joint 
symposium with ASHS under the chairmanship of J. A. Rigney with 
papers on “Plot Technique” by G. F. Potter, E. P. Lana, E. J. Koch, 
W. T. Federer, and W. C. Jacob, and on “Taste Testing and Quality 
Evaluation” by R. A. Bradley and I. J. Koch and D. D. Mason. The 


Wednesday morning session was on “Enumeration Problems in Biology” 
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with H. W. Norton presiding and with papers by C. I. Bliss, C. J. 
Maloney and D. 8. Robson. 


Région Frangaise. A la séance de la Région, le 19 Novembre, 1952, 
4 Paris, le Professeur G. Teissier a parlé sur “Analyse de la variabilité 
chez l’araignée de mer’, et M. M. Lamotte sur “Essai d’emploi de 
quelques méthodes mathematiques en génétique des populations’. 


The Spring meeting of The Biometric Society and The Institute of 
Mathematical Statistics will be held in Washington, D.C., April 29— 
May 1. Please send all contributed papers (talks and abstracts) to 
Boyd Harshbarger, Statistical Laboratory, Virginia Polytechnic Insti- 
tute, Blacksburg, Va., by April 10, 1953. 


In Vol. 8, No. 3 of Biometrics, the publication of the eleven pages of tables in the article, ‘‘Tables for 
Convenient Calculation of Median-Effective Dose (LDso or EDso) and Instructions in their Use’ by 
Carrol S. Weil was made possible by a grant from the Mellon Institute of Industrial Research, Pitts- 
burgh, Pennsylvania. 
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Survival, 89, 90, 271 
Tables, 37, 51, 72, 249, 280, 285, 306, 
353 
Taxonomy, 172, 177 
Termites, 172 
Test of significance, 11, 33, 89, 95, 96, 
112, 170, 222, 275, 295, 384 and 
see non-parametric, sign test 
two-tailed, 36 
Threshold, 388, 389 
Toxicology, 13 
Transformations, 2, 249 
t test, 5, 19, 29, 222, 384 
Variance, 5, 27, 43, 59, 112, 156, 278, 
287, 296, 325, 331, 335, 346, 375 
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